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ABSTRACT 

We present a new iV-body model for the tidal disruption of the Sagittarius (Sgr) dwarf that is capable 
of simultaneously satisfying the majority of angular position, distance, and radial velocity constraints 
imposed by current wide-field surveys of its dynamically young (< 3 Gyr) tidal debris streams. In 
particular, this model resolves the conflicting angular position and radial velocity constraints on 
the Sgr leading tidal stream that have been highlighted in recent years. While the model does not 
reproduce the apparent bifurcation observed in the leading debris stream, recent observational data 
suggest that this bifurcation may represent a constraint on the internal properties of the Sgr dwarf 
rather than the details of its orbit. The key element in the success of this model is the introduction of a 
non-axisymmetric component to the Galactic gravitational potential which can be described in terms 
of a triaxial dark matter halo whose minor/major axis ratio (c/a)$ = 0.72 and intermediate/major 
axis ratio (6/a)$ = 0.99 at radii 20 < r < 60 kpc. The minor/intermediate/major axes of this 
halo lie along the directions (l,b) — (7°,0°), (0°,90°), and (97°, 0°) respectively, corresponding to a 
nearly-oblate ellipsoid whose minor axis is contained within the Galactic disk plane. This particular 
disk/halo orientation is difficult to reconcile within the general context of galactic dynamics (and 
CDM models in particular) , suggesting either that the orientation may have evolved significantly with 
time or that inclusion of other non-axisymmetric components (such as the gravitational influence of 
the Magellanic Clouds) in the model may obviate the need for triaxiality in the dark matter halo. 
The apparent proper motion of Sgr in this model is estimated to be (/Lt/cos&, fib) = (— 2.16, 1.73) mas 
yr _1 , corresponding to a Galactocentric space velocity (U, V, W) = (230, -35, 195) km s _1 . Based on 
the velocity dispersion in the stellar tidal streams, we estimate that Sgr has a current bound mass 
A'^sgr = 2.5^10 x 1O 8 M0. We demonstrate that with simple assumptions about the star formation 
history of Sgr, tidal stripping models naturally give rise to gradients in the metallicity distribution 
function (MDF) along the stellar debris streams similar to those observed in recent studies. These 
models predict a strong evolution in the MDF of the model Sgr dwarf with time, indicating that 
the chemical abundances of stars in Sgr at the present day may be significantly different than the 
abundances of those already contributed to the Galactic stellar halo. We conclude by using the new 
iV-body model to reevaluate previous claims of the association of miscellaneous halo substructure with 
the Sgr dwarf. 

Subject headings: galaxies: individual (Sagittarius) — Galaxy: structure — Galaxy: kinematics and 
dynamics — dark matter 



1. INTRODUCTION 

With the advent of progressively deeper photometric 
surveys, it has become apparent that galactic haloes are 
threaded with the phase-mixed detritus of multiple gen- 
erations of dwarf satellites that have been destroyed by 
the inexorable tides of their host's gravitational poten- 
tial. Surveys have indicated the presence of such streams 
in both the Milky Way (e.g., Majewski et al. 1996; Hclmi 
et al. 1999; Newberg et al. 2002; Rocha-Pinto et al. 
2004; Belokurov et al. 2006; Duffau et al. 2006; Grill- 
mair & Dionatos 2006; Starkenburg et al. 2009) and 
other nearby galaxies (e.g., Shang et al. 1998; Ibata et 
al. 2001a; McConnachie et al. 2003; Kalirai et al. 2006; 
Martinez-Delgado et al. 2008, 2009), possibly constitut- 
ing the primary source of the Galactic stellar halo (e.g., 
Searle & Zinn 1978; Majewski 1993, 2004; Majewski et 
al. 1996; Bullock et al. 2001b; Bullock & Johnston 2005; 
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Font et al. 2006; Bell et al. 2008). 

Perhaps the most prominent such stream is that from 
the Sagittarius dwarf spheroidal galaxy (hereafter Sgr 
dSph) , whose lengthy tidal streams wrap entirely around 
the Milky Way. The first observations of Sgr were pre- 
sented by Ibata et al. (1994), since which time signifi- 
cant observational effort has been invested in detecting 
and characterizing the stellar streams emanating from 
the dwarf. Observations of high-latitude carbon stars 
(Totten & Irwin 1998) and various small-field surveys 
(e.g., Dinescu et al. 2002, and references therein) fu- 
eled a host of early efforts to model the Sgr — Milky 
Way system by Johnston et al. (1995, 1999), Velazquez 
& White (1995), Edelsohn & Elmegreen (1997), Ibata et 
al. (1997), Gomez-Flechoso et al. (1999), Helmi & White 
(2001), Martmez-Delgado et al. (2004). 

In the past seven years however, our understanding of 
the scope and significance of the streams has been revo- 
lutionized by the deep, wide-field views provided by the 
Two Micron All-Sky Survey (2MASS) and Sloan Digi- 
tal Sky Survey (SDSS). In particular, the 2MASS survey 
showed a large population of young, relatively metal-rich 
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Sgr M-giants wrapping 360° or more across the sky (Ma- 
jewski et al. 2003). SDSS observations have shown (Be- 
lokurov et al. 2006) that the debris streamer leading Sgr 
along its orbit continues to be well-defined through the 
North Galactic Pole as it passes over the solar neighbor- 
hood towards the Galactic anticenter. Follow-up spec- 
troscopy (Majewski et al. 2004; Law et al. 2004; Yanny 
et al. 2009) has confirmed that the streams are similarly 
well-defined in radial velocity, and indicate a significant 
evolution in the metallicity distribution function (MDF) 
of the component stellar populations with increasing sep- 
aration from the Sgr dwarf core along the tidal streams 
(e.g., Bellazini et al. 2006; Chou et al. 2007, 2009; 
Monaco et al. 2007; Starkenburg et al. 2009; Carlin 
et al. in prep.). 

In the inner Galaxy, dynamical tracers such as tidal 
streams can provide strong constraints on the mass dis- 
tribution of the baryonic components of the Milky Way 
(e.g., Font et al. 2001; Ibata et al. 2003; Koposov et 
al. 2009). At the larger radii r ~ 15 — 60 kpc traced 
by the orbit of the Sgr dwarf however the gravitational 
potential is dominated by dark matter, and observations 
of the luminous tidal streams can be used to constrain 
the shape, orientation, and mass of the Milky Way's dark 
halo (e.g., Ibata et al. 2001b; Helmi 2004; Johnston et al. 
2005; Law et al. 2005 [hereafter LJM05]; Fellhauer et al. 
2006; Martmez-Delgado et al. 2007). The breadth and 
quality of the recent observational data is such that it 
has highlighted a 'halo conundrum' (variously discussed 
by LJM05; Fellhauer et al. 2006; Martmez-Delgado et 
al. 2007; Newberg et al. 2007; and Yanny et al. 2009): 
In a static axisymmetric Milky Way dark halo no one 
model has been capable of simultaneously reproducing 
both the angular position and distances/radial velocities 
of tidal debris in the Sgr leading arm. Depending on 
which criterion a given study chose to focus, claims have 
been made in favor of an oblate halo (e.g., Johnston et al. 
2005; Martmez-Delgado et al. 2007), an approximately 
spherical halo (e.g., Ibata et al. 2001b; Fellhauer et al. 
2006), and a prolate halo (e.g., Helmi 2004). 

As we demonstrated using massless test-particle orbits 
in Law, Majewski, & Johnston (2009; hereafter LMJ09) 
however, it is possible to resolve this conundrum by 
adopting a triaxial halo in which the minor axis is ap- 
proximately aligned with the line of sight to the Galac- 
tic Center. In this contribution we extend the results 
of LMJ09 by performing detailed ./V-body simulations 
of the tidal disruption of the Sgr dwarf within a range 
of Galactic halo parameterizations that were favored by 
test-particle orbits. In £j3]we outline our basic numerical 
approach, describing the generalized form of the Milky 
Way — Sgr system along with an overview of the obser- 
vational constraints on the Sgr tidal streams. We next 
describe our solution to the problem of constraining the 
basic orbital parameters of the Sgr stream (for a given 
Milky Way model) in <J31 using these various simulations 
to discriminate between different realizations of the un- 
derlying Galactic potential in Sj4] The final parameters 
describing our best-fit model of the Milky Way — Sgr 
system are summarized in Sj5] 

In SjHlwe discuss the best- fit model for the dark matter 
halo of the Milky Way that was derived in <2J and present 
possible theoretical challenges to, alternatives for, and 
methods of testing the triaxial halo model. We demon- 



strate in $7] that by adopting simple assumptions for the 
star formation history of Sgr it is possible to reproduce 
the trends in the MDF observed along the stellar tidal 
streams. We compare the best-fitting model of the Sgr 
streams to other observations of halo substructure that 
have previously been postulated to be associated with 
the Sgr stream in ^ and conclude in <JJ] with a sum- 
mary of the successes and failures of our new model of 
the Milky Way — Sgr system and an overview of the key 
issues that need to be addressed in future work. 

Throughout this paper we adopt the heliocentric Sgr 
coordinate system (Aq,Bq) defined by Majewski et al. 
(2003) in which the longitudinal coordinate A© = 0° 
in the direction of Sgr and increases along trailing tidal 
debris, and £?© is positive towards the orbital pole 
(Z,6) P oic« (274°, -14°). 

2. BASIC APPROACH 
2.1. Numerical Framework 

We adopt a basic formalism similar to that described 
by LJM05 and LMJ09. The Milky Way is described 
by a smooth fixed gravitational potential consisting of 
a Miyamoto-Nagai (1975) disk, Hernquist spheroid, and 
a logarithmic halo. The gravitational potential is there- 
fore specified directly through the isovelocity ellipsoids 
both for consistency with previous studies and because 
these are more analytically accessible than working with 
an underlying mass distribution from which the gravita- 
tional acceleration must be numerically derived. 

The respective contribution of these components is 
given by: 
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where R/r are cylindrical/spherical radii respectively, 
and the various constants C\, C2, C3 are given by 
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As discussed by LMJ09, this form for the dark halo 
potential describes an ellipsoid rotated by an angle <fi 
about the Galactic Z axis, in which q\ and qi are the 
axial fiattenings along the equatorial axes and q z is the 
axial flattening perpendicular to the Galactic disk. We 
adopt the convention that = 0° corresponds to q\/q2 
coincident with the Galactic X/Y axes respectively, and 
increases in the direction of positive Galactic longitude. 
Since only the ratios between the various q are important 
(as increasing/decreasing q\jqilq z in tandem is largely 
degenerate with the total halo mass normalization i/haio), 
we fix (72 = 1.0. 
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We held fixed any parameter that was found not to 
be significantly constrained by Sgr (i.e., that did not 
significantly change the path of Sgr debris when var- 
ied). Following Johnston et al. (1999) and LJM05, 
we take Af d isk = 1.0 x 10 11 M Q , M sphere = 3.4 x 10 10 
M Q , a = 6.5 kpc, b = 0.26 kpc, and c = 0.7 kpc. In 
LJM05 we explored the possibility of scaling the total 
mass of the Galactic disk through the parameter a and 
found no compelling reason to adopt a value other than 
a = 1.0. We fix the distance to the Galactic Center to be 
Rq — 8 kpc (e.g., Reid 1993; Groenewegen et al. 2008); 
in L JM05 we explored the possibility of varying Rq , and 
found that for cases where the Sun lies approximately in 
the plane of the orbit (as is the case for Sgr) any reason- 
able choice of Rq does not have a particularly noticeable 
effect on the quality-of-fit of our simulations. Addition- 
ally, Figure 4 of LJM05 suggests that the ideal choice of 
the scale length of the Galactic halo (rhaio) is reasonably 
insensitive to the axial flattening adopted, and that a 
single fixed value of n^io = 12 kpc should suffice. The 
normalization of the dark halo mass via the scale pa- 
rameter t>haio is therefore entirely specified (for a given 
choice of q\ 1 q2, and 4>) by the requirement that the speed 
of the Local Standard of Rest be «lsr = 220 km s _1 . 
As for Rq, i»lsr is at present relatively unconstrained 
by observations of the Sgr stream (LJM05; although see 
discussion by Majewski et al. 2006; Carlin et al. in prep) 
and we simply fix vlsr to this conventional value. We 
note that if we were instead to adopt the proper mo- 
tion of Sgr A* measured by Reid & Brunthaler (2004) it 
would imply vlsr = 235.6 km s _1 for our adopted value 
of Rq = 8.0 kpc. 

In a recent contribution (LMJ09) we integrated mass- 
less test-particles representing the Sgr dwarf along or- 
bits in the Galactic potential defined by Eqns. [I] -[3] for 
a rang^l of parameters 4> — — 180°, qi — 0.6 — 1.8, 
q z = 0.6 — 1.8. While the accuracy of such simple orbit 
models is necessarily limited because actual tidal debris 
from a massive satellite will deviate slightly from these 
orbits, LMJ09 were able to identify a region of parameter 
space around 4> ~ 90°, q\ « 1.5, q z « 1.25 generally fa- 
vorable to producing a model satellite capable of match- 
ing the observational constraints enumerated in %2.2{ We 
note that while such test-particle orbits technically iden- 
tify a second similarly acceptable region of parameter 
space around 4> ~ 0°, qi « 0.67, q z « 0.83, the physical 
axis ratios of this region are degenerate with the first: 
Both have the short axis approximately aligned with the 
Galactic X axis, and the longest with the Galactic Y. 
We choose to focus on the q z > 1 branch of the solution 
to avoid introducing degeneracy into our simulations. 

Motivated by the successes of LMJ09, we perform N- 
body simulations in Galactic models with a range of pa- 
rameters <j> = 75° - 115°, qi = 1.1 - 1.6, q z = 1.1 - 1.5. 
In each case, Sgr is integrated ~ 8 Gyr backwards in 
time along a simple test-particle orbit, and a 10 5 -particle 
model for the dwarf is integrated forward again in time 
from this point to the present day. The particles in our 

4 Note that LMJ09 only included the results from the range 
qi = 1.0 — 1.8, q z = 1.0—1.8 in their Figures 3 and 4 since much of 
the qi < 1.0, q z < 1.0 parameter space is degenerate with q\ > 1.0, 
q z > 1.0. The non-degenerate parameter space with the minor axis 
oriented along q z was strongly disfavored in the LMJ09 simulations. 



model of Sgr are initially distributed to generate a Plum- 
mer (1911) model 



Sgr 



Sgr,0 



(7) 



where Ms gri o is the initial mass of Sgr and tq — 
(-Msgr,o/10 9 M Q ) 1 / 3 kpc is its scale length. We have as- 
sumed spherical symmetry and isotropy of the orbits 
within the initial Sgr dwarf since: 1) There is no evi- 
dence yet of significant rotation in the main body of the 
satellite, and 2) debris resulting from isotropic orbits will 
cover all phase-space regions accessible to non-isotropic 
satellites. The simulation particles are not associated 
with specific dark or light components, but instead trace 
the total mass of the satellite. The mutual interactions 
between these self-gravitating particles^ in the external 
Milky Way potential are calculated using a self-consistent 
field code (Hernquist & Ostriker 1992). 

2.2. Overview of Observational Constraints 

We summarize below the observational constraints 
that we adopt on the physical properties of the Milky 
Way — Sgr system, each of which is discussed in detail in 
S}3] - HI Generally, we have adopted these constraints be- 
cause we consider them to be the strongest, most unam- 
biguous available. Note that we opt not to constrain the 
model to match the apparent bifurcation in the Sgr lead- 
ing arm described by Belokurov et al. (2006) since recent 
observational data presented by Yanny et al. (2009) sug- 
gest that feature may arise due to substructure within the 
Sgr dwarf rather than the shape of the Galactic gravita- 
tional potential (see discussion in £)4.2I and 15. 2j) . 

While individual constraints are discussed in detail in 
Sj3j we summarize below 

1. Sgr must lie at the angular coordinates, distance, 
and radial velocity described in £13.11 

2. Sgr must have a velocity vector that lies within the 
orbital plane defined by recent trailing tidal debris. 

3. The radial velocity dispersion among trailing tidal 
debris must match that observed by Monaco et al. 
(2007) using high resolution echelle spectroscopy. 

4. The run of radial velocities with orbital longitude 
along trailing tidal debris must match data pre- 
sented by Majewski et al. (2004). 

5. The run of radial velocities with orbital longitude 
along leading tidal debris must match data pre- 
sented by Law et al. (2004; see also LJM05). 

6. The angular location of leading tidal debris must 
match the high surface brightness stream observed 
by Belokurov et al. (2006) in the SDSS. 

7. The angular width of leading tidal debris must 
match that observed by Belokurov et al. (2006). 

5 The gravitational influence of the bound satellite is applied 
to all iV-body particles regardless of whether or not they are still 
bound to the dwarf. 
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It turns out to be easiest to constrain the properties 
of the Sgr dwarf itself first O, since many of them are 
largely independent of the exact form of the Galactic po- 
tential (to within the range explored in this paper), and 
the rest may be constrained contingent upon a specific 
potential. We then use these simulations to discriminate 
between different realizations of the underlying Galactic 
potential in fj4] We note that the only distance infor- 
mation explicitly incorporated in our fitting method is 
that used in setting the distance to the Sgr dwarf itself. 
Our analysis therefore relies only upon the most well- 
constrained observations of the tidal debris (namely, the 
angular coordinates and radial velocities) and is insen- 
sitive to photometric distance calibration uncertainties 
which have hindered confident interpretation of previous 
studies. While we do not fit for the distances to Sgr tidal 
debris explicitly, we shall see in Sj5]and £|7.3l that the final 
model is nonetheless an excellent match to the distances 
derived for the observational data using the most reason- 
able assumptions about distance calibrations. 

3. CONSTRAINING THE PROPERTIES OF THE SGR DSPH 

3.1. Core Properties 

We assume that the Sgr dwarf is located at Galactic co- 
ordinates (l,b) = (5.6°, -14.2°) (Majewski et al. 2003). 
Estimates of the distance to Sgr typically range from 
24 - 28 kpc (see Table 2 of Kunder k Chaboyer 2009 for 
a summary); we adopt the distance scale -Dsgr = 28 kpc 
from Siegel et al. (2007) based on Hubble Space Tele- 
scope (HST)/ ACS analysis of the stellar populations in 
the Sgr core. In LJM05 (which adopted £>s gr = 24 kpc) 
we explored the effect of varying -Dsgr arl d found it to be 
largely degenerate with other parameters (such as Rq) 
and not particularly well-constrained, although models 
generally favored higher values of Z?s gr /i?0. We therefore 
simply fix -Dsgr = 28 kpc to minimize the number of free 
parameters in our model. This places the Sgr dwarf at 
Galactic coordinate^ (X,Y,Z) = (19.0,2.7,-6.9) kpc. 
The radial velocity of Sgr is taken to be v la d = 171 
km s _1 (Ibata et al. 1997), and the dwarf is assumed 
to be presently moving toward the Galactic plane (Irwin 
et al. 1996; Dinescu et al. 2005). 

3.2. Constraining the Orbital Plane 

With the advantages of an internally homogeneous se- 
lection function, M giants drawn from the 2MASS survey 
clearly define the orbital plane of the Sgr stream, espe- 
cially those high latitude M giants free of thick disk con- 
tamination and spanning ~ 150° across the South Galac- 
tic Hemisphere (Majewski et al. 2003). This trailing 
stream is visible as it wraps from the Sgr core (A Q = 0°) 
through the South Galactic Hemisphere towards the 
Galactic anticenter (A Q ~ 150°) — near the position 
of the previous apocenter of the Sgr orbit (A Q ~ 170°) 
. The stream appears to dwindle rapidly at A Q > 140° 
however, largely because the increasing age and decreas- 
ing metallicity (e.g., Chou et al. 2007, 2009) of this older 
tidal debris results in a decline in the prevalence of the 
M-giant tracer population. 

6 We adopt a right-handed Galactic Cartesian coordinate system 
with origin at the Galactic Center. All velocities are given in the 
Galactic Standard of Rest (GSR) frame, with respect to which the 
Sun has a peculiar motion (U, V, W) = (9, 12 + 220, 7) km s" 1 . 



The trend of radial velocities along the trailing stream 
is well-defined (Majewski et al. 2004), enabling us to 
use kinematical information to prune Galactic interlop- 
ers and select a sample of extremely high-likelihood Sgr 
trailing arm stars. We draw such a sample from the Ma- 
jewski et al. (2004) catalog, including all stars in the 
range A© = 0° — 140° for which radial velocity data has 
been obtained^ We fit a fourth-order polynomial to the 
radial velocity data as a function of orbital longitude, 
using an iterative 2.5a rejection algorithm to produce 
our final clean sample of trailing arm stars (Figure [TJ 
left panel) Note that we have not applied an explicit 
distance selection cut to the sample of stars plotted in 
Figure [T] if we were to do so it would remove a majority 
of the crosses from the left-hand panel as these stars lie 
significantly closer than the Sgr trailing arm stream (see, 
e.g., Figure 2 of Majewski et al. 2004). 

The visible and well-defined portion of the trailing arm 
seen in the Southern Galactic Hemisphere 2MASS M gi- 
ants is particularly important because it is dynamically 
young (torn off in the last ^0-2 pericentric passages, 
depending on the assumed mass of Sgr) and has there- 
fore experienced very little differential precession from 
the main body of Sgr (e.g., Helmi et al. 2004; John- 
ston et al. 2005). Observations (Majewski et al. 2003) 
indicate that the trailing debris stream is well-confined 
to a single plane, and numerical simulations show that 
iV-body debris in the range A Q = 0° — 140° have angu- 
lar momentum vectors aligned with that of the present 
satellite to within 1 degree for a wide range of models of 
the Galactic gravitational potential (see, e.g., Figures 4 
and 5 of Johnston et al. 2005). 

The lack of orbital precession in the trailing stream al- 
lows us to determine robustly the instantaneous angular 
momentum vector of Sgr. Using the sample of trailing 
arm stars defined above, we use a GC3 cell count tech- 
nique (Johnston et al. 1996) to define the plane that best 
fits all of the trailing arm stars and the present position 
of the Sgr dwarf. For any given orbital pole (l p , b p ) the 
distance Di of star i from the plane defined by the pole 
is 

Di = h- (fi - r Sgr ) (8) 

where n = (cos 6 P cos / p , cos b v sinl p , sin& p ) is the orbital 
pole vector, r, is the vector from the Galactic Center to 
star i, and rs gr is the vector from the Galactic Center to 
Sgr. 

The best fitting plane is that which minimizes the sum 
D 2 = TiiD 2 over the 108 data points in the clean sample 
of trailing arm stars shown in Figure [T] In Figure [5] we 
plot D as a function of (Z p ,& p ), and demonstrate that 
the best-fitting orbital pole is (Z p , b p ) = (273.8°, -14.5°). 
Unsurprisingly (since both the sample of stars and the 
technique were very similar), this is nearly identical 
to the trailing arm pole calculated by Johnston et al. 
(2005), although it differs slightly from the pole adopted 
by LJM05, who used the pole derived by Majewski et al. 
(2003) from fitting to both leading and trailing Sgr tidal 

7 The sampling criteria of Majewski et al. (2004) required target 
stars to have infrared color J—K s = 1.0—1.1, and be located within 
5 kpc of the Sgr plane defined by Majewski et al. (2003). 

8 We provide an online-only table of the coordi- 
nates and radial velocities in this clean sample at 
|http: //www. astro. Virginia. edu/~srm4n/Sgr| 
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Fig. 1. — Radial velocities for candidate Sgr stream stars in the leading (right-hand panel) and trailing (left-hand panel) tidal streams 
from Law et al. (2004, 2005) and Majewski et al. (2004) respectively as a function of Sgr orbital longitude Aq. Open squares (crosses) 
represent stars meeting (not meeting) our iterative 2.5<r rejection algorithm, the solid lines represent the converged fourth-order polynomial 
fits to the velocity trends. 
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Fig. 2. — GC3 orbital pole plot showing the total value of D for 
a grid of orbital poles. Contours are logarithmic in D and indicate 
the best-fitting orbital pole lies at (Z p ,t> p ) = (273.8° , -14.5°). 

debris. 

3.3. Constraining the Orbital Speed of Sgr 

Since the radial velocity vector w ra( j and orbital pole 
f q3.2j) of Sgr are known, the 3D space velocity of Sgr 
(#s g r = «rad + ^tan, where by definition w rad • u tan = 0) 
can be uniquely specified by the magnitude i>t a n of the 
velocity of Sgr tangential to the line of sight: 



WSgr 



Urad ' n + t/tan ■ ft = 



(9) 



In particular, for each choice of the Galactic halo pa- 
rameters 4>,qi,q z there is a value of v tan that produces 
a trailing tidal stream whose velocities best match those 
shown in Figure [TJ We converge upon this value in a 
manner similar to LJM05 by minimizing the x 2 statistic: 



1 



Xv, trail 



trail 



E 



(WmodolW - Vobs,trail[«])' 



(10) 



where v m odci[z] is determined via linear interpolation of 
the model to the A Q [i] of each of the n^traii = 108 ve- 
locity data points u bs, trail [i]< 

There are two components governing the apparent ve- 
locity width of the stream; physical dispersion and obser- 
vational uncertainty in the derived radial velocities. Ma- 
jewski et al. (2004) find a net dispersion of 12 km s _1 , 
composed of 6 km s _1 intrinsic uncertainty in the veloc- 
ity measurements and an underlying 10 km s _1 physical 
dispersion in the stream. For the purposes of Eqn. [10] it 
is irrelevant what the source of the dispersion is (since 
both will have identical effects governing the offset of 
stars from the orbital path) and we take a v — 12 km s" 1 . 
Rather than perform computationally intensive TV-body 
simulations to constrain t> ta n for every choice of <p, qi,q z , 
we instead use test-particle orbits (see, e.g., LMJ09) to 
estimate the necessary value. Based on ~ 150 TV-body 
simulations with a wide range of Sgr masses and forms 
for the Galactic potential, we find a typical correction 
factor where the best-fit ftan for TV-body models is sys- 
tematically 4 km s _1 slower than that implied by x 2 min- 
imization using such test-particle orbits. 

Considering the example case where tj> = 90° and 
qi = 1.0, in Figure [3] (bottom panel) we illustrate the 
dependence of x 2 trail on w tan for three different values of 
q z . Note that there is a well-defined best value (which we 
call fbest) °f ^tan for each q z \ we estimate the location of 
the minimum by fitting a simple third-order polynomial 
to Xv trail as a function of «tan- The resulting values of 
Ubcst are well-behaved as a function of q z , as illustrated 
in the top panel of Figure |3] 

For each choice of </> and q\ it is possible to construct 
a similar relation. For <f> = 90° the relation between 
«bost and q z is practically unchanged (to within ^1 — 2 
km s^ 1 ) for any choice of qi = 0.8—1.6. For other choices 
of <j) and q\ the deviation of the relation from that shown 
in Figure [3] grows as </> departs from 90°, but for the range 
of (j> = 75 — 115° considered here these relations differ by 
less than 10 km s _1 . 
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FlG. 3. — Bottom panel: Reduced xt trail statistic for the fit of a 
tcst-particlc orbit to trailing arm velocity data as a function of the Sgr 
velocity "Utan tangential to the line of sight for fixed (p — 90°, qi — 1.0 
and three choices of q z — 0.80/1.10/1.60. Solid/dottcd/dashcd lines 
respectively indicate third-order polynomial fits to the discrete data 
points. Top panel: Best-fit Sgr tangential velocity v^ cst (defined by 
minimizing Xv trail) as a function of q z for <fi — 90°, q± — 1.0. The 
solid lino represents a polynomial fit to the discrete measurements. 

3.4. Constraining the Mass of Sgr 

Data in the longitude range A Q = 25° — 90° consist of 
a relatively orderly sample of trailing arm stars as they 
pass through the south Galactic cap moving predomi- 
nantly transverse to our line of sight, and stars in this 
range can provide the most accurate measurement of the 
intrinsic velocity dispersion of the strearrQ- Majewski et 
al. (2004) performed such an analysis, concluding that 
for a reasonable estimate of the observational uncertainty 
that the intrinsic velocity dispersion of the stream is 
0t,,obs — 10.0 ± 1.7 km s" 1 . More recent, higher spectral 
resolution observation^ by Monaco et al. (2007) have 
revised this value slightly to ov bs = 8.3 ± 0.9 km s _1 . 

Since the velocity dispersion in the tidal debris streams 
reflects the mass present within the tidal radius of Sgr 
on the orbit immediately prior to that debris becom- 
ing unbound (see, e.g., Johnston et al. 1996; Johnston 
1998; Peharrubia et al. 2008), it can therefore be used to 
constrain the mass of the dSph (e.g., LJM05). We per- 
formed a series of ~ 160 iV-body simulations in which 
satellites with 10 different initial masses ranging from 
M Sgr ,o = 6.67 x 10 7 - 2.57 x 1O 9 M were integrated 
along reasonable orbits for various lengths of time in dif- 
ferent realizations of the Milky Way gravitational poten- 
tial. The resulting relation between the final bound Sgr 
mass Msgr and the velocity dispersion a v is shown in Fig- 
ure 01 Extrapolating from the best-fit polynomial to the 

9 There is a slight thickening in the distribution of radial veloc- 
ities at A© > 100° because of the angular superposition of tidal 
debris from different epochs (i.e., stripped from Sgr on different 
pericentric passages) that occurs when this debris 'piles up' near 
the previous apogalactic point along the orbit. 

10 Note that while we use the Monaco et al. (2007) velocity 
dispersion to derive the mass of Sgr, we used the Majewski et al. 
(2004) velocity dispersion in Equation llOl bccause this was the value 
appropriate to that more extensive sample of stars used to define 
the velocity trend of the trailing arm. 
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FlG. 4. — Velocity dispersion er v along the Sgr trailing arm (filled 
boxes) as a function of the final bound Sgr mass Afs gr for ~ 160 W-body 
simulations with different realizations of the Galactic potential, Sgr 
initial mass Afggr.Os and interaction time. The hatched area and dashed 
line indicate the observational range cr v — 8.3 ± 0.9 km s _1 measured 
from high-resolution spectra by Monaco et al. (2007). The solid curve 
represents a polynomial fit to the simulated data. 

curve, we conclude that a present bound mass Mg gr = 
x 10 s M© produces the best fit to the observed 
velocity dispersion. Adopting a typical mass-loss history 
(generally similar along reasonable orbits in the range of 
Galactic potentials explored here), this implies that the 
original mass of Sgr was M Sgr ,o = x 10 s Af© if it 

has been orbiting in the Galactic potential for ~ 8 Gyr. 
We adopt Ms gr ,o = 6.4 x 10 8 Af Q , corresponding to an 
initial Sgr scale length ro = 0.85 kpc (see Equation [7]) . 

4. CONSTRAINING THE PROPERTIES OF THE MILKY 
WAY HALO 

We have now constrained the properties of Sgr itself 
as functions of the Galactic parameters cf>,qi,q z . In the 
following subsections we go on to explore within which of 
these potentials the observational characteristics of the 
tidal streams are best reproduced. 

4.1. Sgr Trailing Arm Properties 

As defined in H3.31 we use the statistic \v, trail to quan- 
tify the fit of a given ./V-body simulation to the trailing 
arm velocity trend. While Xv. trail has already been mini- 
mized for a given choice of (j>, q\ , q z by determining the or- 
bital velocity of the satellite, the precise minimum value 
can vary slightly between different halo realizations and 
we therefore also use it as a constraint on the properties 
of the Milky Way. 

4.2. Sgr Leading Arm Properties 

Although both the leading and trailing tidal tails are 
traced by the M-giant population, the location of Sgr 
slightly past pericenter along its orbit in combination 
with foreshortening of the orbital rosettes caused by our 
location 8 kpc from the center of the Milky Way means 
that while the trailing arm extends for ~ 150° across the 
southern Galactic hemisphere the leading arm is clearly 
visible only for ~ 75° (a significant section of which is 
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obscured as it passes behind the Galactic disk/bulge). 
The exact track of leading Sgr debris as it passes through 
the north Galactic cap towards the anticenter is therefore 
uncertain in the 2MASS data since the M-giant tracer 
population is no longer as prominent in this older section 
of the stream. In contrast, SDSS imaging is sensitive to 
older Sgr populations, and clearly shows (Belokurov et 
al. 2006; Newberg et al. 2007; Cole et al. 2008) the Sgr 
leading arm arcing through the northern Galactic cap 
towards re-entry into the Galactic disk in the direction 
of the anticenter. 

Belokurov et al. (2006) initially suggested that there 
may be a bifurcation in the angular position of the 
Sgr stream in the northern Galactic hemisphere, with 
a high surface brightness 'A' stream representing the 
primary wrap of the leading tidal arm passing through 
the coordinates (a, S) ~ (160°, 20°) and a lower surface- 
brightness 'B' stream representing a secondary wrap of 
the trailing tidal arm offset to higher declination around 
{a, S) ~ (160°, 30°) (see also Fellhauer et al. 2006). We 
adopt the positions of the 'A'-stream survey fields of Be- 
lokurov et al. (2006) along this high-surface brightness 
branch as our constraint on the orbital path of leading 
Sgr tidal debris. These 17 pointings are tabulated in 
Table [TJ Note that Belokurov et al. (2006) adopted a 
photometric distance scale calibrated against a distance 
to Sgr of 25 kpc. The distances listed in Table [1] have 
been rescaled to our adopted £*sgr — 28 kpc by multiply- 
ing by a factor of 1.12. 

We note that the dynamical origin of the lower surface- 
brightness B stream is uncertain: Recent observational 
data presented by Yanny et al. (2009) indicate that 
the A and B streams have indistinguishable trends in 
distance, radial velocity, metallicity, and mix of stellar 
types. This suggests that the A and B branches repre- 
sent debris stripped from Sgr at similar times, and that 
the bifurcation may likely be within a single phase of 
tidal debris, originating in substructure within the ini- 
tial Sgr dwarf (see discussion in $5.2\i . Since we do not 
attempt to constrain the internal structure of the Sgr 
dwarf (we focus here on exploring the effects of triax- 
iality in the gravitational potential of the Milky Way, 
and such an effort is beyond the scope of this contribu- 
tion), we therefore do not explicitly constrain our model 
to match the B branch of the bifurcation. 

Numerically, we define the statistic 



xi 



i 



(^modelW - £>obsH) Z 



(11) 



describing the accuracy with which the declination of the 
iV-body model traces the A stream of Belokurov et al. 
(2006); and 



(w<5,modcl — Ws^bsY 



(12) 



describing the accuracy with which the angular width 
of the TV-body model matches the SDSS data at a right 
ascension a — 180° (the non-spherical gravitational po- 
tential gives rise to angular broadening of the leading 
debris stream in this region due to differential preces- 
sion). Based on the data presented by Belokurov et al. 
(2006), we estimate that the la angular width of the 



TABLE 1 
Leading Arm SDSS Constraints 81 



RA 


Dec 


Distance* 3 


(degrees) 


(degrees) 


(kpc) 



215 
210 
205 
200 
195 
190 
185 
180 
175 
170 
165 
150 
1 15 
140 
135 
130 
125 



4 
4 
4 
4 
7 

9.5 
11.2 

13 
13.75 

15 

16 
18.4 

19 
19.4 
19.5 
19.5 
19.5 



51.9 
51.0 
50.0 
47.3 
43.7 
41.1 
38.2 
35.7 
32.5 
30.3 

22.6 " 
21.9 
20.5 



+ 14.6 
2.2 
15.5 
2.8 
+6.0 
4.5 
+ 3.9 
2.7 
+2.9 
1.7 
+5.5 
4.2 
+3.9 
2.5 
+2.7 
1.7 
+3.6 
2.7 
+2.5 
1.1 



+2.3 
1.8 
+3.5 
2.6 
+ 1.8 
1.2 



18.8 



+3.0 
2.3 



a Survey fields are from Belokurov et al. (2006). 

Distances and uncertainties arc taken from Belokurov ct al. (2006) 
and scaled to reflect our choice of distance modulus to the Sgr core. 

stream as seen in SDSS is u^obs = 5.2°, with an un- 
certainty a w = 1.3°. In contrast to Equation [TU] (which 
used A a [i] as the variable coordinate), <5 m odei[*] is deter- 
mined via linear interpolation of the model to the right 
ascension a[i] of each of the ns = 17 SDSS survey fields 
along the main branch of the Sgr stream, each of which 
is assumed to have an angular uncertainty^ a$ — 1.9°. 

While the trend of leading tidal debris is difficult to 
trace with 2MASS imaging data alone, the stream can 
be recovered using spectroscopy to identify the radial ve- 
locity signature of the stream. In Figure [T] we show our 
selection function for stars in the Sgr leading arm with 
K s > 9 (see Ej7.3[) from the observational data presented 
by Law et al. (2004; see also LJM05), using a method 
similar to that adopted for trailing arm stars. This ve- 
locity trend has recently been confirmed by spectroscopic 
observation of a larger number of stars drawn from the 
SDSS imaging of the leading tidal stream (Yanny et 
al. 2009), in contrast with the preliminary conjecture 
by Newberg et al. (2007) that the velocity trend may 
be artificially produced by confusion in the leading arm 
M-giant sample. This gives our final constraint on the 
stream 



Xv.lead 



1 



^"U.lead 



E 



Ddcl U I ^obs.lcad 



erf. 



(13) 



where w m odoi[*] is determined via linear interpolation of 
the model to the orbital longitude A Q [z] of each of the 
?Viead = 94 observational data points in the 2MASS 
leading arm radial velocity sample. As for the trailing 
arm velocity constraint, we adopt an uncertainty a v = 12 
km s _1 for each radial velocity measurement representing 

11 This is a rough estimate of the uncertainty in the centroid of 
the stream at a given right ascension based on the apparent width 
of the stream. 
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a combination of observational uncertainty and intrinsic 
stream width. 

4.3. Combined Constraints 

In general, the parameters <j),q\,q z are not separa- 
ble and must be constrained simultaneously by Xv, trail, 
Xv,icad, XSt and Xw- We combine these four constraints 
into a single statistic: 

X 2 = Xv.trail + Xv.load + x\ + Xw (14) 

and seek the combination of <f), qi , q z that gives the global 
minimum for this combined statistic. 

We performed a series of ~ 500 TV-body simulations 
covering the range 4> = 75° — 115°, q\ = 1.1 — 1.6, q z 
I , I — 1.5 with a grid spacing of 5° in (j>, 0.05 in q z , and 
0.1 in q\. The resulting values of x as a function of 
0, q±, q z are shown in Figure[5]and demonstrate a general 
preference for values of q\ q z in the range <fi rs 90° — 
105°. Interpolating along the grid surface, the optimal fit 
(X = 3.4) is achieved for $ = 97°, q x = 1.38, q z = 1.36. 
We use these parameters for our final best-fit iV-body 
simulation. 

Since this optimal fit has % > 1 there are obviously 
some residual imperfections in the model because dif- 
ferent observational constraints favor different regions of 
parameter space. The greatest mismatch lies in the lead- 
ing arm radial velocities (x-w.lead = 2.7), followed by the 
leading arm angular coordinates (xs = 1-6), the trailing 
arm radial velocities (x«, trail = 1-3), and lastly the an- 
gular width of the leading stream (x w — 0.9). For com- 
parison, the spherical halo model discussed by LJM05 
is clearly a poor fit to the observational data (see Fig- 
ure 10 of LJM05) and has x ~ 9. In general we deem 
simulations with x > 5 to be visibly poor fits. 

The optimal shape for the Galactic dark halo therefore 
appears to be an approximately oblate ellipsoid whose 
short axis is located in the Galactic disk along the direc- 
tion (I, b) = (7°,0°). We note that this differs somewhat 
from that derived by LMJ09, who found a more strongly 
triaxial ellipsoid with = 90°, qi — 1.50, q z = 1.25 by 
matching observed tidal debris to the orbital path of a 
point mass orbiting in the Galactic potential. It is well 
known, however, that actual leading/trailing arm debris 
falls inside/outside of the orbit respectively (e.g., Eyre 
& Binney 2009), meaning that a slightly different shape 
will be derived for the Galactic halo when attempting 
to match simulated debris to observed debris (as done 
in this contribution), rather than a simulated orbit to 
observed debris (as in LMJ09). 

5. MODEL OVERVIEW 

5.1. Model Summary 

In LJM05, we presented possible models for the tidal 
disruption of Sgr under three different assumptions; that 
the Galactic halo was axisymmetric in the disk plane 
and either oblate, spherical, or prolate in the direction 
perpendicular to the disk. Neither of these three mod- 
els were entirely satisfactory as they could not simul- 
taneously reproduce the angular position, distance, and 
radial velocity trends of the leading arm. In contrast, 
our new modelQ orbiting in a non- axisymmetric dark 

12 A complete data file for the dynamical and chemical charac- 
teristics of Sgr debris from this best-fit model is provided on the 



matter halo, reproduces all of these constraints to rea- 
sonable accuracy as demonstrated by Figure [H Indeed, 
we note that even though we didn't explicitly fit to the 
distances of the SDSS leading arm data from Belokurov 
et al. (2006), the new model nonetheless matches them 
to within the uncertainty of the photometric calibration. 
While it is not possible to test the agreement of this new 
model with the distances of the M-giant sample directly 
(due to the evolution of the MDF along the Sgr tidal 
streams), we demonstrate in %7.3\ that the model satis- 
factorily reproduces their apparent magnitude trend. 

Since the Sgr streams frequently overlap themselves 
in all of the observable coordinates, it is useful to iso- 
late specific wraps of the leading/trailing streams in or- 
der to discuss individual wraps more easily. By sorting 
the streams according to angular momentum (relative to 
the Sgr core) and time at which a given debris parti- 
cle became unbound from the model satellite this is eas- 
ily accomplished in the TV-body simulations. We color- 
code debris particles according to the time they became 
unbound from the satellite: Black points are currently 
bound to Sgr, orange points were stripped from the dwarf 
during the last two perigalactic passages (0 - 1.3 Gyr 
ago), magenta points during the previous two (1.3 - 3.2 
Gyr ago), etc. (see Figure [7] for a detailed accounting 
of the color-coding scheme). As a convenient shorthand 
notation, we describe specific wraps of the streams as 
follows (see summary in Figure [8j : 

LI:: Primary wrap of the leading arm of Sgr (i.e., 0° — 
360° angular separation from the dwarf). Traced 
by orange/magenta points near the Sgr dwarf, ma- 
genta/cyan/green points at larger angular separa- 
tions. 

L2:: Secondary wrap of the leading arm of Sgr (i.e., 
360° — 720° angular separation from the dwarf). 
Traced almost entirely by cyan/green points. 

Tl:: Primary wrap of the trailing arm of Sgr (i.e., 0° — 
360° angular separation from the dwarf). Traced 
by orange/magenta points near the Sgr dwarf, ma- 
genta/cyan/green points at larger angular separa- 
tions. 

T2:: Secondary wrap of the trailing arm of Sgr (i.e., 
360° — 720° angular separation from the dwarf). 
Traced almost entirely by cyan/green points. 

The major observational distinction of the present 
model from the oblate-halo model discussed by LJM05 is 
that the primary LI wrap of the leading Sgr tidal stream 
is no longer predicted to pass closer than ~ 10 kpc to the 
Sun; the center of the LI streamer intersects the Galac- 
tic disk at a distance of 13.5 kpc in the direction of the 
Galactic anticenter (see, e.g., Figures [6] and [8]) . This is 
consistent with the absence of net vertical flow observed 
in the solar neighborhood in recent radial velocities sur- 
veys (Seabroke et al. 2008). 

Due to the triaxial Galactic potential, the orbital pe- 
riod and pericenter/apocenter of Sgr varies from orbit 
to orbit (Figure [7]); the mean values of these parame- 
ters are 0.93 Gyr and 19/59 kpc respectively over the 

Web at |http://www. astro,virginia,edu/~srm4n/Sgr 
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Fig. 5. — The combined quality-of-fit criterion x ls plotted as a function of the dark halo axial scales q\ and q z for nine values of the 
rotation angle in the range <j> = 75° — 115°. The vertical stretch varies between panels, but the color map is identical. The best fit is 
achieved for dark green/black regions around q± S3 q z in the central panels <f> = 95° and 100°. The range of \ shown on these plots belie 
the large range in quality of fit of the simulations: While X $ 4 (bright green — black regions) represents a reasonable fit to all of the 
observational criteria, X ^ 10 (blue and magenta regions) are extremely discrepant with at least one criterion. 



past 8 Gyr. The present space velocity of Sgr in this 
model is (U,V,W) = (230, -35, 195) km s" 1 . For our 
assumptions about the distance to Sgr, the Local Stan- 
dard of Rest, and the solar peculiar motion this trans- 
lates to a proper motion of (/x;cos6, = (—2.16,1.73) 
mas yr -1 ([/i Q cos(5, [1$] = [—2.45,-1.30] mas yr -1 in ce- 
lestial coordinates) in the heliocentric rest frame. This 
agrees within 1.9<7 with the proper motion measure- 
ments of Dinescu et al. (2005) who found (/iicosb, ^b) = 
(-2.35 ± 0.20, 2.07 ± 0.20) mas yr" 1 based on South- 
ern Proper Motion Catalog values for 2MASS-selected 
giant star candidates, and within 2. 2c of the values 
(/XiCOS&,Mb) = (-2.62 ± 0.22, 1.87 ±0.19) mas yr" 1 de- 
termined by Pryor et al. (2010). Our tangential velocity 
for Sgr (wtan = 251 km s" 1 ) is also in good agreement 
with the measures of Irwin et al. (1996; Vt an = 250 ± 90 
km s" 1 ) and Ibata et al. (2001b; w tan = 280 ± 20 
km s _1 ). The mass of the Milky Way within 50 kpc 
in this model is 4.4 x 10 11 M Q (i.e., slightly smaller than 
that obtained by Sakamoto et al. 2003), corresponding 
to a halo scaling parameter Whaio = 121.9 km s _1 (see 
Equation [3]) . 



The present mass of the Sgr dwarf in this model is 
-ft^Sgr = 2.5^1^ x 1O 8 M ; taking the apparent magnitude 
of the bound core of Sgr to be my — 3.63 (Majewski et al. 
2003) gives an absolute magnitude My = —13.64 and lu- 
minosity Lsgr = 2.4 x 10 7 L Q for our adopted Sgr distance 
modulus ([m-M] = 17.27; Siegel et al. 2007) implying 
a mass/light ratio of Ms gr /Lsgr = IOI^Mq/L©. This is 
at the low end of the range Ms gr /Lsgr = 14 — 36A/ Q /L0 
estimated by LJM05; this difference reflects both our 
greater adopted distance to Sgr and our choice to use the 
higher-resolution Monaco et al. (2007) measurement of 
the dispersion of the Sgr stream instead of the previous 
Majewski et al. (2004) value. We note, however, that 
the derived mass is somewhat dependent on the struc- 
tural model adopted for the Sgr dwarf. At present, the 
total interaction time of Sgr (and hence the initial mass 
of the satellite) is relatively unconstrained: With rea- 
sonable assumptions for its initial mass the length of the 
tidal streams observed in 2MASS and SDSS suggest that 
Sgr must have been disrupting for at least the last 3 
Gyr (roughly corresponding to orange/magenta-colored 
points in Figure [6]). We have opted to run our simulation 
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FlG. 6. — The best fit TV-body debris model (colored points) in a triaxial halo, overplottcd with observational constraints from 2MASS + SDSS 
(open symbols). The left/right columns plot leading/trailing arm debris respectively. Colors represent the time at which a given debris particle 
became unbound from Sgr; for clarity we plot only debris from the most recent ~ 3 Gyr of tidal disruption (see color-coding key in Figure [71 for 
details) which best correspond to the observational data shown. Open triangles represent SDSS A-strcam data, open squares leading/trailing arm 
2MASS M-giants culled from Figurc[2] The solid 1 inc represents the orbital path of the Sgr dSph core 1.2 Gyr ahead and behind its current position. 
The Sgr dwarf is currently located at orbital longitude/latitude (Aq , Bq) — (0°,0°). 

for a significantly longer time (~ 8 Gyr) in order to make 
concrete predictions about where Sgr debris from earlier 
epochs may be expected; observational confirmation of 
the presence or absence of such older debris streams will 
provide a strong constraint on the actual length of the in- 
teraction time. Future models that more accurately treat 
the dark/light matter separation and are constrained to 
reproduce the bar- like core structure of Sgr (e.g., Lokas 
et al. in prep.) may provide more stringent constraints 
on the Mj L ratio of the dwarf. 

5.2. The Nature of the Bifurcation 

As indicated by Figure [SJ our simulations match the 
high surface-brightness (hereafter HSB) track of the lead- 
ing arm in the 'Field of Streams' region described by 
Belokurov et al. (2006) but do not produce a bifurca- 
tion in the form of a low surface-brightness (hereafter 
LSB) branch offset to higher declination. Fellhauer et al. 
(2006) previously explained this bifurcation as the juxta- 
position of the young (< 4 Gyr) LI wrap of the leading 
arm with the old (> 7.4 Gyr) T2 wrap of the trailing 
armJ3 This explanation is not viable for the current 

13 Note that Fellhauer et al. (2006) call the equivalent of their 
L1/L2/T1/T2 streams the A/C/D/B streams respectively. 




Time (Gyr) 

FlG. 7. — Distance of the model Sgr dwarf from the Galactic 
center over its ~ 8 Gyr orbital history. Time is in Gyr before 
the present day. Different color-codings represent different orbital 
phases on which debris was lost in the Af-body model (i.e., orange 
debris on ./V-body plots throughout this paper correspond to debris 
lost on the orange colored orbits shown here, etc.). 
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FlG. 8. — Diagram showing the four wraps of the model Sgr stream in Galactocentric Cartesian coordinates (XqCi ^Gc)- The Ll/Tl wraps 
represent the primary leading/trailing tidal streams respectively, while the L2/T2 wraps represent the secondary leading/trailing tidal streams of 
debris which has been wrapped > 360° around the Milky Way from Sgr. Colors of individual debris particles represent the time at which the 
particle became unbound from the Sgr dwarf (sec Figure[7]for a key). The solid black line represents the orbital path of the Sgr dwarf core, filled 
black circles are drawn every 0.2 Gyr along the orbit. The X indicates the location of the Galactic Center, while the circled X indicates the location 
of the Sun. 



model: As illustrated in Figure HI the distance/radial 
velocity trends of the LI (black points) and T2 (green 
points) streams are markedly different, in contrast to 
Yanny et al. (2009) who demonstrate that the distance 
and radial velocity trends of the LSB and HSB streams 
are indistinguishable. 

However, we note that the Yanny et al. (2009) obser- 
vations are not reproduced well by the Fellhauer et al. 
(2006) model either: This model predicts a factor ~ 2 
difference in distance between the LSB and HSB streams 
at a = 120° that is not observed. Additionally, given the 
strong evolution in the MDF with increasing separation 
A© from the Sgr dwarf (e.g., Chou et al. 2007) we should 
not expect the metallicity of the LSB and HSB streams 
to be similar (as observed by Yanny et al. 2009) if they 
represent extremely different dynamical ages of tidal de- 
bris. Certainly, if the LSB stream were significantly older 
than the HSB we should not expect the LSB stream to 
contain as many M-giants as the HSB stream since M- 
giants are a young tracer population. 

A-body models of Sgr tidal destruction also predict 
that the young Tl trailing arm should pass through the 
Field of Streams region. In both the model presented 
here and the model described by Fellhauer et al. (2006), 
the Tl wrap traverses this region from a ~ 120° to 
a ~ 220° with a distance varying from ~ 40 kpc to ~ 20 



kpc, and radial velocity varying from ~ 125 km s to 
~ —125 km s _1 . Such tidal debris (grey points in Fig- 
ure [9]) is not observed in the Field of Streams, suggesting 
that the Sgr tidal streams are insufficiently long to pro- 
duce Tl trailing tidal debris in this region of sky (or at 
least, that the stellar density in the stream is declining 
sufficiently rapidly that it has not been detected in the 
SDSS). This places a limit on the length of the Sgr tidal 
stream similar to that derived from other observational 
constraints that the Sgr tails represent < 3 Gyr worth of 
tidal debris. If the Sgr tails are so young, we should not 
expect any debris corresponding to earlier wraps (i.e., 
T2/L2) to exist and it is correspondingly impossible to 
interpret the bifurcation as the juxtaposition of Sgr tidal 
debris torn from Sgr at significantly different times. 

We therefore favor the explanation proffered by Yanny 
et al. (2009) that both the LSB and HSB streams may 
represent debris stripped from Sgr at similar times. In- 
deed, we note that models of Sgr orbiting in a non- 
spherical Galactic potential exhibit precessional broad- 
ening of the angular width of the leading debris stream 
in this region which manifests as a low surface bright- 
ness 'spray' to higher declinations similar to that of 
the LSB stream (this effect is increasingly noticeable for 
haloes with stronger departures from spherical symme- 
try). This presents an alternative scenario wherein the 
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FlG. 9. — Predicted locations, distances, and radial velocities of Sgr tidal debris in the 'Field of Streams' region described by Bclokurov et al. 
(2006). This figure is shown in celestial (a, S) coordinates instead of the more natural (Aq,Sq) to aid comparison to the figures of Bclokurov et 
al. (2006) and Fcllhaucr et al. (2006). Note that the color-coding in this figure is different than that adopted in all other TV-body debris plots 
(e.g., Figure [8j: Black/dark grey points represent Ll/Tl leading/trailing arm debris, blue/dark green points represent L2/T2 leading/trailing arm 
debris. 

ratiofS (c/a)<s ~ 0.72, (b/a)$ - 0.78 (Hayashi et al. 
2007). While there is a tendency for haloes to be more 
spherical at larger radii (Allgood et al. 2006; Hayashi 
et al. 2007), the halo of Milky Way- like systems at 
a distance of ~ 50 kpc is expected to be prolate with 
X$ « 0.88 (based on the Via Lactea simulation of Dic- 
mand et al. 2007; Kuhlen et al. 2007). 

Obscrvationally, typical haloes appear to be rounder 
than predicted by collisionless CDM simulations, as in- 
dicated by gravitational lensing mass models (e.g., Hoek- 
straetal. 2004; Mandelbaum et al. 2006; Evans & Bridle 
2009), polar ring galaxies (e.g., Iodice et al. 2003), and 
X-ray isophotal studies (e.g., Buote et al. 2002; although 
c.f. Diehl & Statler 2007; Flores et al. 2007). This 
possible mismatch is perhaps unsurprising because the 
growth of baryonic structures within these haloes tends 
to make them considerably rounder (e.g., Dcbattista et 
al. 2008; Valluri et al. 2009), inflating the principal axis 
ratios by as much as ~ 0.2 — 0.4 within the virial radius 
(Kazantzidis et al. 2004). However, it is difficult to con- 
strain the actual triaxiality of individual haloes because 
galaxies are projected onto the two-dimensional plane of 
the sky and observations are relatively insensitive to the 
mass distribution along the line of sight (although c.f. 

14 We adopt a convention in which axial ratios such as c/a are 
given with a subscript of <E> (p) to indicate that they are measured 
with respect to the potential (density) contours. 



bifurcation may arise due to substructure within the ini- 
tial Sgr dwarf, such as a smaller dwarf satellite that may 
have been bound to Sgr when it fell into the gravitational 
potential of the Milky Way (see, e.g., Chou et al. 2009; 
Law et al. in prep) , or simply to anisotropy of the stellar 
orbits within the original dwarf which might be expected 
to give rise to tidal debris occupying only a subset of the 
angular width covered by our isotropic iV-body model. 

The detailed nature of the low surface brightness 
branch in the SDSS view of the Sgr stream is there- 
fore unknown at present, and is not possible to reconcile 
with iV-body models that fit the unambiguous primary 
leading and trailing tidal streams illustrated in Figure [51 
While this feature may ultimately provide a constraint 
on the internal structure of the initial dSph satellite, ex- 
ploration of such structure is beyond the scope of this 
contribution. 



6. DISCUSSION: THE SHAPE AND ORIENTATION OF THE 
MILKY WAY DARK MATTER HALO 

Cold dark matter (CDM) models of galaxy formation 
generally predict that dark haloes should be mildly triax- 
ial (e.g., Bullock 2002; Jing & Suto 2002; Bailin & Stein- 
metz 2005; Allgood et al. 2006; and references therein), 

l—(b 2 /a 2 ) 

predominantly prolate systems (T$ = i-^/a*)® ~ 0.66; 
Allgood et al. 2006) with characteristic central axis 
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Corless et al. 2009). 

It is therefore noteworthy that for the only galaxy 
where we have a fully three-dimensional view from the 
interior of the mass distribution (i.e., the Milky Way) it 
appears f t!4.3[) that the best-fitting axial ratio at radii 
20 < r < 60 kpc describes an almost perfectly oblate 
ellipsoid with minor/major axis ratio (c/a)$ = 0.72, in- 
termediate/major axis ratio (6/a)$ = 0.99, triaxiality 
T<j, = 0.06, and the minor axis located in the plane of 
the Galactic disk in the direction (I, b) — (7°, 0°). In this 
section we discuss the physical interpretation of these re- 
sults, along with possible effects that could influence our 
derivation of the halo shape. 

6.1. Physical Validity of the Triaxial Halo Model 
6.1.1. Tracing the Underlying Mass Distribution 

In i 32.ll we introduced flattening directly into the con- 
tours of the gravitational potential because this was an- 
alytically tractable and consistent with the method of 
previous studies (in particular LJM05). It is well-known 
however that as the ellipticity of the potential becomes 
more extreme, the corresponding mass distribution re- 
quired to generate it may become unphysical (i.e., have 
regions of negative mass). In order for any model of the 
Galactic halo to be believed, it must correspond to a 
physically realizable density distribution. 

In Figure [TU] we show the corresponding contours of 
the mass distribution underlying our best-fit model for 
the Milky Way (as calculated by taking the Laplacian 
of the gravitational potential). As expected, the mass 
distribution is less spherical than the isopotential sur- 
face; the contours of the dark matter halo are approx- 
imately elliptical at small radii (i.e., within ~ 40 kpc) 
with (c/a) p = 0.44, (b/a) p = 0.97, and T p = 0.07. Al- 
though some pinching of the dark matter density con- 
tours is apparent at larger radii the mass distribution is 
everywhere positive. Figure [10] also demonstrates that 
while addition of the baryonic mass component signifi- 
cantly changes the net density distribution at small radii 
and near the Galactic disk plane, the density distribution 
at larger radii probed by the orbit of Sgr is dominated 
by the dark halo component. 

We additionally confirmed the validity of our results by 
constructing a triaxial NFW (Navarro et al. 1996) halo 
model in which the axial flattcnings were introduced in a 
more physically motivated manner directly into the mass 
distribution and where the resulting accelerations were 
computed using the approximate form of the correspond- 
ing potential given by Lee & Suto (2003). As expected, 
there are slight differences in the path of satellites or- 
biting in this NFW halo and in our best-fit logarithmic 
potential, but these differences are minor and likely due 
to the fact that we did not fine tune the NFW model. 

6.1.2. Orientation of Principal Axes 

While the mass distribution underlying our best-fit 
model of the Galactic halo is physically reasonable, the 
orientation of the halo (i.e., strongly non-axisymmetric 
in the plane of the Galactic disk, with the angular mo- 
mentum vector of the disk aligned with the intermedi- 
ate axis of the halo) is unexpected because orbits about 
intermediate axes tend to be unstable (e.g., Adams et 
al. 2007). This is not a concern for stellar orbits in 



the Galactic disk because the gravitational potential at 
small radii r < 20 kpc is dominated by the (nearly) ax- 
isymmetric baryonic mass distribution rather than the 
triaxial dark halo. Indeed, preliminary results of orbit 
integration in the gravitational potential derived here in- 
dicate that loop orbits are permitted, suggesting that a 
self-consistent disk may still be supported (R. Johnson 
and K.V. Johnston, private communication, 2010). How- 
ever, this situation begs the question of why the Galactic 
disk formed in its present orientation unless the princi- 
pal axes of the dark halo either twist substantially with 
radius (in contrast to the general prediction of Hayashi 
et al. 2007) or evolve with time. While we cannot rule 
out the possibility that rotation of the dark halo has con- 
tributed in some regard to our results, we neglect such 
effects because typical haloes are expected to rotate ex- 
tremely slowly (e.g., Bullock et al. 2001a; Herbert- Fort 
et al. 2007, 2008; Bett et al. 2009). 

Typically, simulations of disk galaxy formation suggest 
(e.g., Debattista et al. 2008) that the minor axes of the 
dark halo and the baryonic disk should be approximately 
aligned (to within an average of ~ 30°; see discussion 
by Warren et al. 1992; van den Bosch 2002; Bett et al. 
2009) because both are formed from accreted matter with 
similar specific angular momenta. If a disk were to form 
in an elliptical halo potential such as that proposed in 
£14.31 simulations suggest that the disk too should become 
elongated and non-axisymmetric (see, e.g., Gerhard & 
Vietri 1986; Debattista et al. 2008). It will therefore be 
intriguing to investigate whether the relative orientation 
of the host dark matter halo and baryonic disk proposed 
here can be strongly motivated within the current CDM 
paradigm. 

It is worth asking as well whether such an inclined 
oblate halo can explain the observed distribution of 
Galactic satellites, which simulations suggest (Zentner 
et al. 2005) should be distributed along the major axis 
of the host halo (i.e., have a pole aligned with the minor 
axis). The Fornax-Leo- Sculptor great circle has a pole 
(I, b) = (135°, -2.9°) (Lynden-Bell 1982; see also Kunkel 
& Demers 1976; Majewski 1994; Fusi Pecci et al. 1995), 
and a recent update (Metz et al. 2007, 2009) including 
the faint new dwarfs discovered in the SDSS gives a re- 
vised satellite pole (l 7 b) = (159.7,-6.8). The updated 
pole of Metz et al. (2009) lies within ~ 30° of the minor 
axis we derive for the Galactic halo, but it is difficult to 
determine whether this similarity is simply an unrelated 
coincidenceP^I 

6.2. Alternatives to the Triaxial Halo Model 

In LJM05, we conducted an extensive parameter space 
search in an effort to find a combination of Milky Way 
+ Sgr parameters capable of simultaneously matching 
both the observed angular position and distance/radial 
velocity trends of leading tidal debris. In particular, we 
explored the effect of varying the proper motion of Sgr 
(i>tan), the distance to Sgr (Ds g r), the distance from the 
Sun to the Galactic Center (i?©), the total mass of the 
Galactic disk (via the scaling parameter a in Equation 
[T]) , the dark halo scale length (d) , the dark halo flattening 

15 If the probability distribution of the satellite pole were uni- 
form across the entire sky, the chance for it to be coincidentally 
aligned with (l,b) = (7°,0°) to within 30° is ~ 13%. 
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Fig. 10. — Density maps viewed along different axes for the best-fit Galactic halo model in an (X 1 , Y' , Z) coordinate system aligned with 
the principal axes of the halo. Density maps are shown for the dark halo alone, and for the sum of all mass components (i.e., logarithmic 
halo, Miyamoto-Nagai disk, and Hcrnquist spheroid). The greyscale/red contours are logarithmic and with the same scale in all panels. 
Note that while the axisymmetric density profile of the disk dominates in the disk plane (Z = 0), at only 1 kpc above the plane the halo 
contribution produces noticeably elliptical contours in the X' — Y' plane at radii r > 20 kpc. The dashed green circles indicate the mean 
apocenter/pericenter of the Sgr orbit. 



perpendicular to the disk (q z ), and the total mass scale 
of the Milky Way as parameterized by the rotation speed 
of the Local Standard of Rest (i>lsr)- 

Although the Galactic disk contributes significantly to 
the overall flattening of our adopted Galactic potential 
at small radii (e.g., Figure [TO)) , it has only a minor effect 
at the range of radii (r ~ 20 — 60 kpc) probed by the 
orbit of Sgr. As demonstrated by Figure [TTJ the dark 
halo in our model is the dominant contributor to the 
total gravitational potential at these radii. Indeed, if 
we decrease the mass of the Galactic disk by 50% (by 
setting a = 0.5 in Equation!!) and increase the total halo 
mass correspondingly (so that the Local Standard of Rest 
remains constant at i>lsr — 220 km s^ 1 ), we can repeat 
the parameter space search described in 21 Scaling the 
proper motion i> ta n °f Sgr to again best-fit the trailing- 
arm velocity trend, we find optimal values for qi,q z ,<p 
almost identical to those derived previously (although 
the overall x 2 of the best-fit model is not as good as for 
the case a = 1.0). We are therefore confident that the 
details of the Galactic disk model do not significantly 
bias our results. 

Similarly, the details of the triaxial Galactic bar (which 
we did not include in our model) and stellar halo are un- 
likely to change our results. Neither are sufficiently mas- 
sive to have a great effect on the gravitational potential at 
distances r > 20 kpc, and neither have orientations akin 
to that which we derive for the dark halo. The major axis 
of the bar is thought to lie within ~ 15° - 20° of the X- 



axis (e.g., Blitz & Spergel 1991; Nakada et al. 1991; Mor- 
ris & Serabyn 1996; Babusiaux & Gilmore 2005), similar 
to the minor axis of the dark halo. While the stellar 
halo is significantly triaxial ([c/a] p ps 0.65, [b/a] p = 0.75, 
T p = 0.76; Newberg & Yanny 2006; see also Larsen & 
Humphreys 1996; Xu et al. 2006) with a major axis ly- 
ing approximately in the Galactic plane (~ 50 — 70° from 
the X axis) the minor axis is not apparently contained 
within the X — Y plane but located ^13° from the Z 
axis. 

Throughout the entirety of the L JM05 parameter space 
search, the only parameter that was found to have a sig- 
nificant effect on the angular position / distance / radial ve- 
locity trend of the Sgr leading arm was the flattening of 
the dark halo q z . Although no single choice of q z can si- 
multaneously reproduce all observational constraints, we 
have demonstrated (see also LJM09) that a fully triaxial 
model in which the short axis is approximately aligned 
with the Galactic X axis is capable of doing so. 

Such a triaxial halo is therefore obviously attractive in 
the sense that only iV-body models in such a halo suc- 
ceed in modeling the Sgr - Milky Way system where pre- 
vious efforts (e.g., Helmi et al. 2004; LJM05; Fellhauer 
et al. 2006; Martmez-Delgado et al. 2007) have failed. 
Despite the extensive exploration of parameter space un- 
dertaken by ourselves (this paper; LJM05; LMJ09) and 
other groups (e.g., Helmi et al. 2004; Fellhauer et al. 
2006; Martmez-Delgado et al. 2007) however, the space 
explored to date is far from exhaustive and it is possible 
that the triaxial halo model may prove to be simply a 
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Fig. 11. — Left-hand panel: Contribution of the dark halo, stellar disk, and stellar spheroid to our model Galactic rotation curve as 
a function of radius in the Galactic disk plane (for purposes of illustration the halo is taken to be axisymmetric with gi = <?2 = 1.0). 
The vertical dotted lines denote the mean pericenter/apocenter of the Sgr orbit. Right-hand panel: Acceleration of the Sgr dwarf in the 
gravitational potential of the Milky Way as a function of time over its assumed 8 Gyr orbital history. The dark halo, stellar disk, stellar 
spheroid, and total contributions to the rotation curve and acceleration of Sgr are indicated by dashed, long-dashed, dot-dashed, and solid 
lines respectively in both panels. 



numerical crutch that mimics the effect of some as-yet 
unidentified alternative. In the following sections, we 
briefly discuss a few possibilities that may obviate the 
need for triaxiality to explain the extant data. 

6.2.1. Alternative Formulations for the Gravitational 
Potential 

While we have focused on a logarithmic formulation of 
the Galactic gravitational potential, other forms could 
also reasonably be adopted. In LJM05 we explored ax- 
isymmetric NFW (Navarro et al. 1996) models and were 
unable to resolve the halo conundrum. As demonstrated 
in Pi6.1.1l however. fully triaxial NFW models can give rise 
to orbits for Sgr that match the observed angular preces- 
sion and distance/radial velocity trends of leading tidal 
debris. We find that orbits within logarithmic and NFW 
haloes are sufficiently similar that there is no reason to 
differentiate between the two for our present purposes. 

Fellhauer et al. (2006) claim that a mass distribution 
of the form specified by Dehnen & Binney (1998) can pro- 
duce some effects similar to prolateness in a logarithmic 
halo. However, using such a model these authors were 
nevertheless unable to reproduce the distance trend of 
leading tidal debris as seen in the SDSS (Belokurov et 
al. 2006) and concluded that their Set D of proper mo- 
tions for Sgr (i.e., those derived by LJM05) in a logarith- 
mic halo provided the best overall fit to the observational 
data then available. We therefore find no reason to pre- 
fer a Dehnen & Binney (1998) model over the triaxial 
logarithmic halo model presented here. 

6.2.2. Evolution in the Orbit 

Although our model dSph has been orbiting in the 
Galactic potential for ~ 8 Gyr, we have not considered 
the influence of any evolution in either the orbit of Sgr 
or the underlying gravitational potential of the Milky 
Way. We expect there was actually evolution in both 



the depth, and the shape of the potential (e.g., Kuhlen 
et al. 2007). While we expect the streams to be indepen- 
dent of past evolution in the potential of the Milky Way 
(since tidal streams respond adiabatically to changes in 
their host potential; see discussion by Penarrubia et al. 
2006) some evolution has certainly occurred in the orbit 
of Sgr since it fell into the Milky Way. Dynamical fric- 
tion and other such evolutionary effects therefore reflect 
an uncertainty in our models, although we note that (as 
discussed by LJM05) such effects alone are unlikely to be 
able to resolve the halo conundrum. Further investiga- 
tion of these effects must (among other things) properly 
account for the relative distribution and mass-loss his- 
tory of both the dark and baryonic matter in Sgr, and is 
therefore beyond the scope of this contribution. 

6.2.3. Dwarf — Tail Gravitational Interaction 

As discussed by Choi et al. (2007), the gravitational 
influence of the bound satellite on debris in the tidal tails 
can shift the trajectories of actual tidal debris away from 
the orbital path traced by the bound core of the dSph. 
For sufficiently massive satellites (greater than ~ 0.1% 
the virial mass of the host), this effect can be apprecia- 
ble and degenerate with that caused by flattening (or tri- 
axiality) in the underlying gravitational potential of the 
Milky Way. However, we note that despite claims to the 
contrary (Choi et al. 2007) our simulations fully account 
for this effect; the self-gravity of the bound satellite is 
applied to all TV-body particles at all times regardless of 
whether or not they are in the Sgr core or in the tidal 
tails. 

6.2.4. Gravitational Influence of the Magellanic Clouds 

The Large Magellanic Cloud (LMC) is ~ 50 kpc dis- 
tant in the direction (l,b) = (280.5°, -32.8°) (van der 
Marel et al. 2002), placing it within ~ 0.5 kpc of the 
Galactic Y — Z plane. Based on great-circle fits to the 
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Magellanic stream, the pole of the system appears to be 
(7 pj & p ) = (188.5°, -7.5°) (Nidever et al. 2008). That 
is, the pole of the Magellanic stream is aligned with the 
short axis of the Galactic dark halo as derived in this con- 
tribution to within ~ 1° in Galactic longitudeF^I This 
suggests two possibilities: (1) The Magellanic Clouds 
may have fallen into the Milky Way along the plane of 
the long and intermediate axes of the dark matter halo 
(i.e., the alignment is a consequence of the dark matter 
distribution), or (2) the Magellanic Clouds may exert a 
significant gravitational acceleration on the Sgr dwarf, 
with the mass of the Clouds at large radii in the Y — Z 
plane mimicking the effect of an extended dark matter 
distribution in this plane (i.e., the alignment is the cause 
of the apparent dark matter distribution). 

A precise description of the second scenario would re- 
quire a comprehensive model for the orbital and mass- 
loss history of the Magellanic Clouds, but it is possible 
to understand the magnitude of their gravitational influ- 
ence on the path of Sgr tidal debris by integrating mass- 
less test-particles along the Sgr orbit in a Milky Way with 
a spherical dark halo and simplified model of the LMC. 
We first test the effect of introducing a fixed point-mass 
into the gravitational potential of the Milky Way at the 
present location of the LMC, whose mass is allowed to 
vary from 1 - 10% of the total mass of the Milky Way 
within 50 kpc (i.e., M LMC = 6 x 1O 9 M - 6 x 10 10 M Q ). 
As illustrated in Figure [12j a fixed LMC with a mass 
10% that of the Milky Way can produce a significant de- 
viation on the order of 20° in the orbital path of leading 
Sgr tidal debris. 

Of course the LMC is not fixed in space, but describes 
an orbit about the Milky Way confined approximately to 
the Y — Z plane. We therefore develop a second, slightly 
more realistic model in which the LMC is described in a 
time-integrated sense by a ring of mass at a radius of 50 
kpc in the Y — Z plane. The gravitational acceleration at 
an arbitrary point in space resulting from such a configu- 
ration is described using an analytical approximation to 
the numerically evaluated elliptical integrals (see deriva- 
tion for the electrostatic case by Zypman 2005). This 
second model produces even more noticeable deviations 
from the original orbital path, with appreciable shifts in 
both the leading arm angular position and radial veloci- 
ties expected for models in the range of masses 1 — 10% 
that of the Milky Way. 

The detailed effect of the Magellanic Clouds on the or- 
bit of Sgr is extremely uncertain. Both our point-mass 
and ring-mass models are obviously oversimplifications, 
neglect the SMC entirely, and assume that the gravita- 
tional influence of the LMC has been fixed over the en- 
tire interaction history of the Sgr dwarf (~ 8 Gyr in the 
current model). Since the Clouds are currently thought 
to be near the pericenter of their orbit, for the major- 
ity of the lifetime of Sgr they likely lay at greater dis- 
tances where their gravitational force would be smaller. 
Indeed, some recent lines of investigation suggest that 
they may have only recently fallen into the gravitational 
potential of the Milky Way for the first time (e.g., Besla 
et al. 2007), drastically limiting the time during which 
they could have torqued the orbit of the Sgr dwarf. Of 

16 They differ by 7.5° in Galactic latitude, but the symmetry 
axes of the dark halo were fixed to lie at b = 0°. 



the various orbital paths shown in Figure I12[ only the 
triaxial halo model actually succeeds in reproducing all 
of the observational data. However, given the myriad 
uncertainties in the Milky Way - LMC - Sgr system the 
important conclusion to draw is that gravitational petur- 
bations from the Magellanic Clouds may be sufficient to 
torque the leading arm of the Sgr stream significantly, 
especially if the effect is enhanced by the response of the 
Milky Way dark matter halo itself to the passage of the 
Clouds (see, e.g., Weinberg 1998). 

6.2.5. Alternative gravity 

Modified Newtonian gravity (MOND; see Milgrom 
1983; Bekenstein 2004) on Galactic scales has been pro- 
posed as an alternative theory to explain the flat rotation 
curves typical of galaxies without recourse to an extended 
halo of dark matter. Read & Moore (2005) demonstrated 
that MOND was an equally viable alternative to the 
oblate dark matter halo model presented by LJM05, ca- 
pable of fitting the angular position of the leading tidal 
debris but not the radial velocities. Now that a dark 
matter model has been able to successfully match both 
observational constraints, it is worthwhile to investigate 
whether a similarly successful MOND model can be con- 
structed. Since the key to resolving the halo conundrum 
appears to have been the inclusion of a strongly non- 
axisymmetric component to the gravitational potential, 
it is not immediately obvious how such a potential could 
be produced in a MOND model by the predominantly 
axisymmetric distribution of baryonic matter typically 
adopted for the Milky Way (although c.f. Wu et al. 2008; 
Widrow 2008). 

6.3. Testing the Model 

Fortunately, the model presented here for the shape, 
size, and orientation of the Galactic gravitational poten- 
tial is easily testable; it must be capable of explaining 
the orbital paths of other Galactic satellites that simi- 
larly provide constraints on the mass distribution in the 
Milky Way. Typically, previous studies of tidal debris 
systems in the Galactic halo have only explored the ef- 
fects of changing the axial scalelength perpendicular to 
the Galactic disk (i.e., q z ), and it is unclear whether sat- 
isfactory solutions for other satellites can be found in 
such a strongly non-axisymmetric potential as that sug- 
gested here. In particular, can realistic orbits be found 
to match such systems as the Monoceros tidal stream 
(e.g., Peharrubia et al. 2005), the Magellanic CloudJH 
(e.g., Besla et al. 2007; Ruzickaet al. 2007), the "orphan 
stream" (Belokurov et al. 2007), the GD-1 stream (Ko- 
posov et al. 2009), the newfound Cetus Polar Stream 
(Newberg et al. 2009), and tidally disrupting Galac- 
tic globular clusters such as Pal 5 (Odenkirchen et al. 
2009)? Similarly, is a triaxial Galactic halo consistent 
with observations of the Galactic HI disk (e.g., Oiling 
& Merrifield 2000; Kalberla et al. 2005; and references 
therein)? It is intriguing that Saha et al. (2009) discuss 

17 While simulations generally favor haloes which are extended 
along the Galactic Z axis (similar to those in the triaxial model), 
preliminary results indicate that the Magellanic Clouds are rela- 
tively insensitive to the distribution of mass along the Galactic Y 
axis and can neither confirm nor deny the triaxial model presented 
here (G. Besla, priv. comm.). 
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FlG. 12. — Effect of the LMC on the angular path, distance, and line of sight velocity of the leading/trailing Sgr orbit (loft- and right-hand panels 
respectively). Observational data (open symbols) arc as in Figure [6] The cyan line represents the path of the Sgr dwarf in a simple spherical halo 
model, green solid/dashed lines the path in a spherical halo with a point-mass LMC that is 1%/10% of the Milky Way mass, and red solid/dashed 
lines the path in a spherical halo with a ring-model LMC that is 1%/10% of the Milky Way mass. It is difficult to sec the solid green line in 
places since the cyan line lies almost on top of it. Note the strong deviations of the models incorporating the LMC from the predictions of a simple 
spherical halo model, especially for larger values of the LMC mass. The black line indicates the path of the dwarf in the triaxial halo model for 
comparison. 



evidence to suggest that the dark matter halo must be 
similarly non-axisymmetric in the Galactic XY plane in 
order to explain the asymmetric flaring of the disk (al- 
though c.f. Narayan et al. 2005). Future observations of 
hypervelocity stars at large distances may also provide 
key diagnostics (e.g., Gnedin et al. 2005; Yu & Madau 
2007; Perets et al. 2009). 

7. EXPLAINING THE EVOLUTION OF THE MDF ALONG 
THE TIDAL STREAMS 

Recently, numerous measurements have been made 
(e.g., Vivas et al. 2005; Chou et al. 2007, 2009; Monaco 
et al. 2007; Yanny et al. 2009; Starkenburg et al. 2009) 
of the metallicity distribution function (MDF) along the 
Sgr stream, although a uniform picture of the evolution 
of the MDF is difficult to construct because systematic 
differences in the metallicity sensitivity between stellar 
populations traced by various surveys can bias the re- 
sulting trends. However, significant evidence is starting 
to emerge for a gradient of the MDF with Sgr orbital 
longitude. Bcllazini et al. (2006) for instance observed 
that the ratio of blue horizontal branch to red clump 
stars changes between the main body of Sgr and the stel- 
lar streams (demonstrating an age/metallicity gradient), 
and similar claims of such evolution have been made by 
Monaco et al. (2007) who found that while the core 



of Sgr has a mean metallicity ([Fe/H]) = -0.35 ± 0.19 
the tidal streams have a mean metallicity of ([Fe/H]) = 
—0.70 ± 0.16. The most comprehensive picture pre- 
sented to date however has been that of Chou et al. 
(2007), whose uniform selection technique demonstrates 
that the MDF of the Sgr M-giant population continu- 
ously evolves from a median [Fe/H] ~ —0.4 in the Sgr 
core to [Fe/H] ~ —1.1 dex over the ~ 180° length of the 
leading tidal arm. 

In this section, we attempt to reproduce the gradient 
in the MDF by making physically motivated assumptions 
about the star formation history of the Sgr dwarf. We 
consider only the M-giant population observed by Chou 
et al. (2007) and Monaco et al. (2007) in order to en- 
sure a uniform selection function at all orbital longitudes. 
While Chou et al. (2009) have recently extended their 
earlier work (Chou et al. 2007) by additionally measur- 
ing the evolution in the abundance patterns of titanium 
(Ti), yttrium (Y), and lanthanum (La), we do not incor- 
porate these additional elements into our analysis. 

7.1. Enrichment Model 

Generally, we expect that star formation in Sgr is likely 
to have occurred near the bottom of its gravitational po- 
tential well, and that stars thus formed in the core would 
gradually diffuse outwards with time. Remnants of such 
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stellar populations would tend to persist in the Sgr core 
until the present day and we therefore assume that the 
stellar populations present in the Sgr streams correspond 
to counterparts in the Sgr core. These Sgr core popu- 
lations have been traced using high-precision HST-ACS 
photometry by Siegel et al. (2007), who find five major 
star formation episodes: 

Population A:: M54 metal-poor population. 
[Fe/H] = —1.7 ±0.1, age 13 ± 1 Gyr. 

Population B:: Sgr metal-poor population. 
[Fe/H] = -1.2 ± 0.1, age 11 ± 1 Gyr. 

Population C:: Sgr intermediate population. 
Three sub-populations: 
CI: [Fe/H] = -0.65 ± 0.1, age 7 ± 0.6 Gyr, 
C2: [Fe/H] = -0.45 ± 0.1, age 6 ± 0.6 Gyr, 
C3: [Fe/H] = -0.3 ± 0.1, age 5 ± 0.6 Gyr. 

Population D:: Sgr young population. 

[Fe/H] = -0.05 ± 0.15, age 2.5 ± 0.5 Gyr. 

Population E:: Sgr very young population. 

[Fe/H] = +0.55 ± 0.05, age 0.75 ± 0.1 Gyr. 

Younger stellar populations might be realistically as- 
sumed to be more tightly bound than older populations, 
which have had more time to diffuse outward in the 
dwarf. Motivated by the presence of such population 
gradients within other Galactic dSphs (e.g., Tolstoy et 
al. 2004) for which the more metal-rich populations are 
more centrally concentrated and have lower velocity dis- 
persion, we tag particles in our best-fit model of the Sgr 
dwarf to correspond to different stellar populations ac- 
cording to their total energy in the pre-interaction dSph. 
The energy of an individual particle in the satellite is 
given by 

m = \vW + m (is) 

where v[i] is the velocity and the internal potential 
energy of particle i relative to the satellite. For conve- 
nience, we define E* to be an order-sorted version of the 
energy distribution, ranging from (most strongly bound 
particle) to 1 (least bound particle) with a uniform dis- 
tribution in the range 0—1. 

We assign particles with a given range of E* to each 
population as illustrated in Figure the locations of 
the energy divisions (and hence the total number of par- 
ticles assigned to each burst of star formation) are free 
parameters. Once a particle has been assigned to a given 
population, we assign it a metallicity and age based on 
a Gaussian random selection about the mean and stan- 
dard deviation of the population given in the list above, 
and randomly choose an infrared color uniformly from 
the range J — K s = 1.0 — 1.1. Adopting Marigo et al. 
(2008) isochronef3 with a Chabrier (2001) IMF for each 
of the simulated populations, we calculate the absolute 
K s magnitude for each simulated particle, adopting a 
50%/50% mix of RGB and AGB stars. This particu- 
lar mix of RGB/ AGB stars is not intended to match 
any particular distribution, but to illustrate the general 

18 We use the interface made available by L. Girardi at 
http: //stev. oapd. inaf . it/ cgi-bin / cmd 
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Fig. 13. — Metallicity assigned to all particles in the simulation as 
a function of their initial energy E* in the original satellite. Dashed 
lines indicate the adopted energy d ivisi ons between the popula- 
tions, as named in the list in Section l7.ll Colors of points indicate 
the pericentric passage on which the particle becomes stripped from 
the satellite (as in Figure |B). 

range of magnitudes expected for stars in each evolu- 
tionary phase. We note that not all of the Siegel et 
al. (2007) populations will produce M-giants: Popula- 
tion A in particular is too old and metal- poorH All this 
means, however, is that we would not expect regions of 
the stellar streams populated by these stars to be traced 
by M-giants of those metallicities, and we therefore do 
not calculate K s magnitudes for these particles. 

Strictly speaking, the enrichment model presented here 
is unphysical since it assumes that all of the stellar popu- 
lations were present in Sgr at all times during its ~ 8 Gyr 
interaction history with the Milky Way, despite the fact 
that some of these populations didn't form until many 
Gyr later. This is an inescapable consequence of the N- 
body method adopted in this paper — particles cannot 
readily be introduced into the simulation at arbitrary 
times. Ideally, a fully hydrodynamical simulation could 
build a gas-phase component into the model dwarf, and 
permit this to form stars at the appropriate times. How 
to retain an appreciable gas supply in Sgr until only a 
Gyr ago is a challenging question in its own right, but 
claims of a detection of gas removed from Sgr as recently 
as ~ 0.4 Gyr ago by its most recent passage through the 
extended Galactic disk have been made by Putman et al. 
(2002). Since the iV-body particles only represent dy- 
namical tracers however, we content ourselves with not- 
ing that no appreciable quantities of any given stellar 
population are tidally stripped from the dwarf in this 
model before they are thought to have formed (see hori- 
zontal dashed lines in Figure [T4|). 

7.2. Matching the Observed MDF 

19 Indeed (as discussed by Siegel et al. 2007), the M54 metal- 
poor population may have formed separately from Sgr and later 
been accreted by the dwarf, since M54 has a distinct stellar pop- 
ulation and radial profile and Sgr is nucleated even ignoring this 
population (e.g., Monaco et al. 2005). 
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Tidal stripping of a satellite on an eccentric orbit such 
as that of Sgr occurs nearly impulsively, with the ma- 
jority of tidal debris stripped in the brief time interval 
that the satellite is passing through the pericenter of its 
orbit. Although the overall trend with successive peri- 
centric passages is to remove stars at successively smaller 
radii or energies (i.e., 'onion-skin' style stripping mod- 
els, see discussion by Martmez-Delgado et al. 2004), 
each mass loss event does not impose a sharp tidal cut- 
off on the stellar populations of the dwarf. Rather, the 
rapidly changing tidal radius, the impulsive nature of 
the stripping, and the dependence of the effectiveness 
of tidal stripping on the orientation of individual stel- 
lar orbits (see, e.g., Read et al. 2006a) means that on 
any given pericentric passage Galactic tides reach deeply 
into the satellite and can unbind stars from a variety of 
populations. This does not completely evacuate the en- 
ergy/radius range of particles in the satellite, and the 
remaining stars rapidly relax back to a Plummer model 
over a dynamical time (~ 0.1 — 0.2 Gyr; see, e.g., dis- 
cussion by Peharrubia et al. 2008, 2009). Despite the 
sharp divisions that we introduced in the distribution of 
stellar populations within the model satellite, there is 
significant overlap in the populations expelled from the 
satellite into the tidal tails over the lifetime of the dwarf. 
This is illustrated in Figure [TJ] While particles assigned 
to Population A are almost entirely stripped from the 
dwarf during the hrst two to three orbital periods, stars 
in Population B are present in many successive epochs 
of tidal debris and still have a residual bound popula- 
tion (black points at i U nbound = —1) in the dwarf at the 
present day. Stars in Population D are present in only 
the most recent tidal debris. The trend of stellar popu- 
lations along the tidal streams is further affected by the 
angular phase smearing caused by the differential orbital 
velocities of stripped stars from each epoch of tidal de- 
bris. 

By adjusting the relative sizes of each of the stel- 
lar populations (i.e., sliding the vertical lines in Fig- 
ures [12] and [T4"|) it is possible to adjust the relative bal- 
ance of stars of a given metallicity present in each de- 
bris era (i.e., point 'color'), and thereby as a function 
of orbital longitude along the Sgr tidal streams. In 
Figure [15] we plot the metallicity of stars observed by 
Chou et al. (2007) and Monaco et al. (2007; their 
SARG subsample) as a function of orbital longitude, 
overplotted with the mean metallicity of the simulated 
Sgr stream using three different assumptions about the 
energy ranges E* assigned to each of the Populations 
A-E. The simplest possible assignment in which Popula- 
tions A/B/C/D/E correspond to uniformly wide energy 
ranges E* = 1.0-0.8/0.8-0.6/0.6-0.4/0.4-0.2/0.2-0.0 
respectively fails to match the observational data. As 
illustrated by the red curves in Figure [15l the highest- 
metallicity debris is over-represented in the regions of 
the stream nearest the satellite (i.e., A Q sw 250° and 
100° in the leading/trailing streams respectively), and 
Population C is over-represented relative to Populations 
A and B at larger angular separations (i.e., the model 
stream is too metal-rich at A Q 30° and 250° in the 
leading/ trailing streams respectively). It is possible to 
suppress the presence of the youngest stellar populations 
in the Sgr stream by reducing the range in E* assigned to 
them, and the green curves in Figure [15] show the result 
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FlG. 14. — Time a given particle is stripped from the model dSph as a 
function of initial energy E* (values of i un b oun d = — 1 indicate particles 
bound to the dSph at the present day). Colors correspond to specific of 
orbits in the lifetime of the dwarf (see Fig. [7] for key). Vertical dashed 
lines indicate the divisions assigned between stellar populations (as in 
Figure [l3t ■ horizontal dashed lines indicate the formation time of each 
population. Note that while mass loss is always occurring, the majority 
(i.e., corresponding to the highest point density) is concentrated into 
relatively brief periods surrounding pericentric passages. 

of adopting energy ranges E* = 1.0 — 0.8/0.8 — 0.6/0.6 — 
0.2/0.2-0.04/0.04-0.0 for Populations A-E respectively. 
While this improves the agreement in the regions of the 
tidal tails nearest Sgr, the mean model metallicity is still 
too high at larger angular separations. 

We find that a reasonable fit to both the low- and 
high-angular separation observational data is achieved 
with a choice of energy ranges E* = 1.0 — 0.8/0.8 — 
0.45/0.45-0.38/0.38-0.29/0.29-0.20/0.20-0.04/0.04- 
0.00 for Populations A/B/C1/C2/C3/D/E respectively 
(black lines in Figure [T5]) . Given the large spread of 
the observational measurements and the large number 
of assumptions made by our debris-tagging method, a 
more precise determination of the best-fit choice of en- 
ergy ranges is not warranted. Rather, we conclude that 
for some choice of the energy ranges it is possible to 
broadly fit the observational trends: The mean metallic- 
ity of the model Sgr stream decreases from [Fe/H]~ —0.4 
in the core to [Fe/H]~ —1.0 at large angular separations. 

The above choice of energy ranges corresponds to a rea- 
sonably strong ~ 2 dex metallicity gradient with radius 
in the initial model Sgr dSph (Figured!)] left-hand panel), 
consistent with the prediction of Chou et al. (2007), 
and slightly stronger than observed in the Sculptor dSph 
by Tolstoy et al. (2004). The strong initial gradient is 
largely erased by the process of tidal stripping, leaving a 
relatively weak metallicity gradient ~ 0.8 dex remaining 
by the present-day (Figure [T51 right-hand panel). The 
MDF of the bound core of the model Sgr dwarf (Figure 
\W\ right-inset black histogram) is therefore significantly 
different than both the MDF of the original dSph (Fig- 
ure [16] left-inset black histogram) and the MDF of the 
stellar tidal debris that Sgr has contributed to the Galac- 
tic halo (Figure [151 right-inset grey histogram) . If such 
strong tidal evolution of the MDF is common for Galactic 
satellites, it may partially explain the observed chemical 
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FlG. 15. — Mctallicity as a function of orbital longitude along the leading/trailing (right/left panels respectively) Sgr tidal streams for the Chou 
ct al. (2007) and Monaco ct al. (2007) observational samples. Different symbols are used to indicate stars in the various leading/trailing arms 
subsamplcs described by these authors. In the middle and lower panels we ovcrplot the mean mctallicity (solid lines) for ./V-body particles in the 
most recent 5 Gyr of tidal debris. The range of typical mctallicitics about this average is indicated by the dashed lines. The red/green lines 
respectively in the middle panels indicate two methods of debris-tagging which do not result in a trend matching the observational data, while the 
black lines in the lower panel generally reproduce the observational trends (see discussion in Q7.21 . 
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Fig. 16. — Metallicity as a function of radius in the pre-interaction (left-hand panel) and present day (right-hand panel) Sgr dSph. Filled 
squares represent the mean value of [Fc/H] at a given radius, the vertical bars represent the lc spread about this mean value at each 
radius. Black points represent [Fe/H] for particles bound to the Sgr dwarf, grey points represent unbound stars in the tidal streams. Inset 
histograms at the left and right edges represent the MDF of the intial and present-day Sgr dwarf respectively: Note the difference between 
the present-day MDF of bound particles (black histogram) and particles in the tidal streams (grey histogram). 
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abundance mismatch between Galactic halo stars and the 
present-day dSph population (e.g., Geisler et al. 2005; 
and references therein). Rather than requiring the early 
accretion onto the Milky Way of metal-poor, high [a/Fe] 
satellites that are completely disrupted by the present 
day (e.g., Font et al. 2006), the metal-poor halo stars 
may simply have originated in the dSphs with which we 
are familiar, and in which few traces of the originally size- 
able metal-poor populations now remain (e.g., Majewski 
et al. 2002; Munoz et al. 2006; Chou et al. 2009). 

7.3. K magnitudes of M giants 

One of the challenges in attempting to trace Sgr de- 
bris in the 2MASS M-giant view originally provided by 
Majewski et al. (2003) was that as the metallicity of 
stars evolves along the stream, so does the zeropoint of 
the photometric distance calibration. Since the distance 
calibration was tied to the relatively metal-rich Sgr core, 
distances become systematically more inaccurate with in- 
creasing separation from the core. Now that we have a 
model for the evolution of metallicity along the streams, 
we can perform the comparison of simulations and ob- 
servational data in a purely observational apparent K s 
magnitude space. 

In Figure [T71 we plot the apparent K s magnitude, an- 
gular position (A Q , 73 Q ), and radial velocity ^gsr for the 
simulated Sgr dwarf, overlaid with the Chou et al. (2007) 
and Monaco et al. (2007) data. We note that the TV-body 
model generally succeeds in reproducing the radial veloc- 
ity trends indicated by these data, although the Chou et 
al. (2007) trailing/leading arm sample of stars wrapped 
into the North/South Galactic Hemisphere (crosses/filled 
squares in Figure ITT]) may constitute evidence for a more 
massive Sgr progenitor in order to produce extended 
wraps of the stream that are broader in radial veloc- 
ity space. The TV-body models also suggest that some of 
the observational data may originate elsewhere than the 
Sgr stream: A few of the crosses, three-five filled circles, 
and two open circles in Figure [17] have radial velocities 
strongly inconsistent both with the other data points and 
the A-body models. 

The agreement between the apparent K s magnitudes 
of the A-body model and the Monaco et al. (2007) data 
appears to be relatively good. Unfortunately it is diffi- 
cult to evaluate this agreement for the Chou et al. (2007) 
data since these authors focussed on the brightest avail- 
able M-giants, which are expected to be biased several 
or more a brighter than the mean for the stream (see 
discussion by Chou et al. 2009). A more instructive 
comparison is to the M-giant survey of Majewski et al. 
(2003). In Figure[l8]we plot the apparent K s magnitude 
for the simulated Sgr dwarf, overlaid with the original 
M-giant sample from Majewski et al. (2003; i.e., similar 
to their Fig. 8) and the subsample with radial velocity 
data which we culled from Law et al. (2004) and Majew- 
ski et al. (2004) in Figure Q] (open boxes in Figure US]). 
The K s magnitudes derived from the model match the 
larger M-giant sample of the dynamically young sections 
of the tidal streams (i.e., A < 150° for the trailing arm, 
and A > 220° for the leading arm) well: Over 95% of 
the radial velocity subsample have K s magnitudes con- 
sistent with the A-body model. The remaining 5% may 
be Galactic interlopers at smaller distances whose col- 
ors and radial velocities happened to be coincident with 



those of the Sgr stream, although a few located around 
{A Q ,K S ) = (280°, 9.5) may be associated with the Tl 
wrap of the trailing arm stream instead of the fainter 
(and more distant) LI stream since the radial velocity 
trends of the LI and Tl wraps are expected to cross in 
this region of the sky. 

There are a few regions of stars in Figure [T5] that are 
not explained by the model however. In particular there 
is a distribution of M-giants at K s ~ 13 — 14 at all or- 
bital longitudes that are not reproduced by the A-body 
model. At A Q ~ 180° - 200° these faint M-giants ar- 
guably appear to join smoothly to the trailing arm trend 
(see also Newberg et al. 2003), suggesting (Pakzad et 
al. 2004) that the trailing arm may not turn around at 
this location to approach closer to the Sun again (as the 
A-body model does) but instead forms a wrap around 
the entire sky at relatively constant K s ~ 13 — 14. Such 
a wrap is strongly inconsistent with the predictions of all 
A-body models to date and, if genuine, may be evidence 
for a strong evolution in the orbital path of Sgr ~ 3 Gyr 
ago. 

7.4. Caveats 
7.4.1. Relation to Other Tracers 

The picture presented here of the evolution in the mean 
metallicity along the Sgr tidal arms is biased towards 
more metal-rich stars, since it has been fitted to obser- 
vations of intrinsically metal-rich M-giant tracers. Other 
tracers such as blue horizontal branch (BHB) stars (e.g., 
Yanny et al. 2009), RR Lyrae (RRL; Vivas et al. 2005, 
Prior et al. 2009), or K-giants (Starkenburg et al. 2009; 
see also Dohm-Palmer et al. 2001) might reasonably 
be expected to differ from the trends calculated here. 
It is unclear, however, how stars more metal-poor than 
[Fe/H] ~ -2 (as detected by Yanny et al. 2009) might 
arise in large numbers in relatively young sections of the 
Sgr stream (i.e., the LI arm at A ~ 250°) since no rem- 
nants of populations so metal poor remain in the core of 
Sgr today (Siegel et al. 2007). If this metal-poor stellar 
population has indeed recently been completely stripped 
from Sgr it may provide a useful constraint upon future 
models of the distribution of stellar populations within 
the initial Sgr dSph. We discuss this issue further in SJH] 

7.4.2. Uniqueness 

The fact that we have been able to find a method for 
tagging particles in the initial satellite to correspond to 
the stellar populations seen in the tidal streams should 
not be taken to mean that this structure of the initial 
satellite is uniquely correct, that the optimal ranges in 
E* found for each star formation episode A-E correspond 
exactly to the total stellar mass produced in the bursts, 
or even that the method is a particularly good descrip- 
tion of the tidal stripping of specific stellar populations 
within a real Galactic dSph. The mass loss history of the 
satellite is a strong function of the adopted model, and 
different values for the total mass, scale length, or inter- 
action time of the dwarf will produce different amounts of 
overlap between various eras of tidal debris, and various 
ranges in orbital energy within the original dwarf. Sim- 
ilarly, if we were to adopt a different physical model for 
Sgr (e.g., a cuspy Hernquist profile, or a more physically 
motivated two-component King-profile stellar core com- 



22 



Evolution of Sgr in a Triaxial Halo 



us 
tU 
ID 
U 
&0 
OJ 
X) 

o 
CD 



40 
20 

■20 
-40 



~i — I — I — i — i — I — I — I — I — i — i — i — I — i — I — I — r 




i — H 


i i i i 


i i i i 


mm" 






MM 












-H- 


1 k M 




1 1 M 




n — i — I — i — i — i — r 




300 200 100 
A<-, (degrees) 



300 200 100 
A (degrees) 



FlG. 17. — TV-body debris particles (colored points) for the best-fit model, overlaid with stars (see key in Figure Il5t for which mctallicity 
information has been acquired by Chou ct al. (2007) and Monaco ct al. (2007). Note that for purposes of comparison we plot more epochs of tidal 
debris here than shown in Figurc[6] For clarity we omit the Sgr core sample of Chou ct al. (2007; i.e., open triangles) from this figure since they 
would be indistinguishable from the black-colored ./V-body debris points. The upper/lower 'branches' of the model stream at slightly different K s 
magnitudes represent the RGB and AGB evolutionary phases respectively. The solid line in the upper and middle panels indicates the orbital path 
of the Sgr dwarf core. 



bined with an extended NFW dark matter halo) , or con- 
sider the possible conversion of an originally disky dwarf 
to a dSph via tidal stirring (e.g., Mayer et al. 2006) the 
mass loss history could also be very different (see, e.g., 
Read et al. 2006ab). The specific distribution and mass 
normalization of individual stellar populations is there- 
fore contingent on the adopted model and should not be 
interpreted as a general result, but can be refined as more 
observational data become available. 

8. DISCUSSION: RELATION TO POSSIBLE DETECTIONS 
OF SGR TIDAL DEBRIS 

The model presented in this paper has the attractive 
quality that it provides a reasonable match to the major 
dynamical characteristics of Sgr tidal debris measured 
from wide-field surveys. It is worthwhile also investigat- 
ing the degree of correspondence between the ./V-body 
model and various narrow-field observations of subgiant 
stars, K-giants, and RRL that have been postulated to 
trace the Sgr tidal stream. 

In Figure [19] we plot the A^-body model of the Sgr 
stream overlaid by a selection of detections from recent 
narrow-field surveys (this comparison is not exhaustive, 
see also discussion by Majewski et al. 2003). We dis- 
cuss each data set in detail below. Note that we do not 



consider possible associations with globular clusters or 
faint dwarf satellites; these are treated separately in a 
forthcoming contribution (Law et al. in prep.). 

8.1. SDSS Stripe 82 

Cole et al. (2008) used a maximum likelihood method 
to characterize the spatial properties of F-turnoff stars in 
stripe 82 of the SDSS. They observed tidal debris located 
at {a,S,R) = (31.37 o ±0.26°,0.0°,29.22±0.20kpc), cor- 
responding to Sgr longitudinal coordinates (A@,Bq) = 
(105.2°, -1. 15°). As illustrated by Figure Q2] (filled cir- 
cles) the Tl wrap of the ./V-body model is an extremely 
good match to these F-turnoff stars in both angular posi- 
tion and distance (matching to within 0.7° in Bq, and 1 
kpc in distance). Cole et al. (2008) measured the width 
of the stream perpendicular to the line of sight at this A Q 
to have a FWHM = 6.74±0.06 kpc. The A^-body stream 
is slightly fatter with FWHM = 7.9 kpc, but given the 
single component mass-follows-light model employed for 
Sgr the agreement with the apparent width of a specific 
stellar population is extremely good. 

The radial velocity of a sample of RRL in Stripe 82 
has recently been measured by Watkins et al. (2009), 
who detect a component with i>gsr = —130 km s _1 , 
which is matched by the mean radial velocity of the Tl 
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FlG. 18. — Apparent K s magnitude as a function of orbital longitude for leading and trailing tidal arms. Colored points represent debris from 
the iV-body model (for purposes of comparison we plot more epochs of tidal debris here than shown in Figure [6}, overlaid with small black points 
representing M-giants from the sample of Majewski et al. (2003; see their Figure 8). We also overlay the subsamplc of M-giants with radial velocities 
(open boxes) pruned in Figure [Tl 

stream at this position to within 5 km s _1 . Watkins et 
al. (2009) also measure the metallicity of these RRL, 
finding that ([Fe/H]) = -1.41 ± 0.19. This value is ap- 
preciably more metal-poor than either the M-giant sam- 
ple or the mean of the TV-body model at this longitude 
([Fe/H] ~ —0.7). This difference in measured metallicity 
emphasizes the strong sensitivity of derived metallicities 
of a stellar stream to the biases imposed by the type of 
stellar tracer used (e.g., metal-poor populations do not 
generally form M giants while metal-rich populations do 
not generally form RRL). Thus, until an unbiased assess- 
ment of the MDF along the Sgr stream is obtained, mod- 
els such as that presented in ^71 which was constrained 
by a single type of metallicity-biased population tracer, 
must be considered as merely a demonstration of how 
MDF-variations can originate. 



8.2. RRL and K-giants at A Q - 270° - 300° 

Vivas et al. (2005) published VLT spectroscopy of 
16 RRL from the QUEST survey in the range A Q ~ 
270° -300°. All but 1-2 of these stars, (which have highly 
discrepant velocities, at vqsr = —106 and —197 km s _1 ) 
are in agreement with the angular position, distance and 
radial velocity trend of the LI arm of the iV-body model 
(see Figure [T9l open boxes). Starkenburg et al. (2009) 



have also detected an overdensity of 5 K-giant stars at 
similar orbital longitude that are well-matched to the 
distance and radial velocity of the iV-body model stream 
(Figure [191 skeletal triangles), although they appear to 
be located towards the edge of the stream in -B©E3 

Both samples of stars are well-reproduced by the ^V- 
body model, although we note that the Vivas et al. 
(2005) and Starkenburg et al. (2009) samples have 
mean metallicity ([Fe/H]) = —1.77 and —1.68 respec- 
tively. These metallicities are entirely consistent with 
each other, but significantly more metal-poor than the 
M-giant sample, which has ([Fe/H]) ~ —0.9. As was 
the case for SDSS Stripe 82, this reinforces the conclu- 
sion that stellar populations significantly more metal- 
poor than the M-giant subsample exist at a range of 
orbital longitudes throughout the Sgr tidal streams. As 
expected on the basis of dynamical age, the RRL of Vivas 
et al. (2005) which correspond to tidal debris ~ 1.5 — 3 
Gyr old (i.e., magenta points in Figure IT9|) are slightly 
more metal-poor than the RRL in SDSS Stripe 82 which 
correspond to ~ — 1.5 Gyr old debris. 

20 Note that the star added to the Starkenburg et al. (2009) 
group using the '4dist' < 0.08 criterion is not as well-matched to 
the TV-body stream as the other 5 stars, with the former discrepant 
by ~ 60 km s _1 . 
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FlG. 19. — Model (colored points) of the Sgr stream overlaid with various detections of potential Sgr tidal debris. We include 2 Gyr more of tidal 
debris in this plot (green points) than in previous plots (e.g., Figure [6j to indicate the possible path of debris if Sgr has been orbiting the Milky 
Way for more than ~ 5 Gyr. Overlaid on the TV-body model is observational data from Vivas ct al. (2005; open boxes), the SDSS 'Stripe-82' (Cole 
ct al. 2008; Watkins et al. 2009; filled circle), Duffau et al. (2006; open circle), 'VOD' observations from Vivas ct al. (2008; filled boxes), K-giant 
stars from Starkcnburg ct al. (2009; skeletal triangles) and Kundu et al. (2002; filled triangles), RR-Lyrae stars from Prior ct al. (2009; open 



triangles), and red clump stars from Majcwski ct al. (1999; crosses). Note that at Aq 



260° 



■ 300° it is difficult to distinguish the open circle of 



Duffau et al. (2006; because it lies atop the open triangles of Prior et al. 2009 in Aq vs distance space) 
et al. (2009; because they lie amidst the open boxes of Vivas et al. 2005). 



and the skeletal triangles of Starkcnburg 



We note that an additional sample of metal-poor K- 
giant stars in a similar longitude range was also pointed 
out by Kundu et al. (2002; filled triangles in Figure [T9]) . 
While the L2 /T2 wraps of the iV-body stream match the 
radial velocity trend of these stars, the angular positions 
are discrepant by up to 20° for L2 and 30° for T2. At a 
typical distance ~ 5 kpc the Kundu et al. (2002) stars are 
also much closer than expected for the T2 arm (d ~ 50 
kpc), but may be consistent with the L2 arm (d ~ 10 
kpc). If these stars genuinely belong to the Sgr stream, 
it is most likely that they belong to the L2 wrap, and may 
suggest that the angular coordinates and/or distance of 
this wrap in the iV-body model require adjustment. 

8.3. Virgo: The VSS and the VOD 

Recent years have witnessed a wealth of substructure 
discovered in the direction of Virgo, which we follow 
Martinez-Delgado et al. (2007) in classifying into the 
Virgo Stellar Stream (VSS) and the Virgo Overdensity 
(VOD). 

The VSS (as traced by QUEST RRL; Duffau et al. 
2006) is a relatively well-defined feature that lies at an 
angular position (Aq,Bq) = (265°, 14°), i.e., projected 
on the outer edge of the L1/L2 streams and significantly 



off the T1/T2 streams in sky position. At a distance of 
19 kpc, it is roughly coincident with the L2/T1 stream, 
but has a radial velocity (ugsr = 99.8 km s _1 ; Duffau et 
al. 2006) that is discrepant with these wraps by 220/94 
km a" 1 respectively. Curiously, the proper motion of the 
VSS ([n a cosS, us] = [-4.85 ± 0.85, 0.28 ±0.85]; Casetti- 
Dinescu et al. 2009) is reproduced by the Tl stream to 
within Ict (\jj, a cos6, us] = [-4.18 ±0.41, 0.04±0.25]), but 
given the angular position and radial velocity mismatch 
we conclude that the VSS is unlikely to be associated 
with the Sgr stream. 

In contrast, the VOD (Juric et al. 2008) is a much 
less well-defined, diffuse clump of stars around (I, b) — 
(300°, 65°) that covers more than ~ 1000 deg 2 on the sky. 
The possible association of the VOD with Sgr has been 
addressed in detail by Martinez-Delgado et al. (2007), 
who noted that it was roughly coincident with tidal de- 
bris from models of the Sgr dwarf (LJM05) disrupting in 
an axisymmetric, oblate Galactic halo. 

In the triaxial halo model derived in this contribution, 
we find that the model Sgr stream arcs significantly over 
the solar neighborhood, and tidal debris in the LI tidal 
stream at the angular coordinates of the VOD lies much 



Law and Majewski 2010 



25 



too far away (~ 46 kpc for the Sgr stream versus ^5 — 17 
kpc for the VOD) to be associated (see Figure H9l filled 
squares). Similarly, both the Tl and T2 wraps lie too 
far away to produce the VOD, in addition to lying in the 
opposite angular direction Bq. 

Interestingly, this angular position and distance does 
correspond closely to the predicted location of the sec- 
ondary L2 wrap of leading Sgr tidal debris, suggesting 
that the VOD may be evidence for old tidal debris torn 
from Sgr > 3 Gyr ago. If such an identification is cor- 
rect, we should expect the radial velocity of the VOD 
to be ^gsr = —131 ± 22 km s _1 . Recent observational 
work on QUEST RRL in the direction of the VOD by 
Vivas et al. (2008) derived radial velocity signatures 
at v GS R = +215, -49, -171 km s" 1 . The first (and by 
far most significant) two of these velocity structures are 
strongly discrepant with the expected signal by ~ 345 
and 80 km s _1 respectively. The third structure is pos- 
sibly consistent with the velocity signature expected of 
the L2 Sgr stream, differing by only ~ 40 km s" 1 in this 
part of the simulated debris where the iV-body stream 
is ill-constrained. Only 7% of the QUEST RRL (i.e., 
3/43) are contained in this «gsr = —171 km s" 1 velocity 
structure however, and we therefore concur with the con- 
clusion of Vivas et al. (2008) that while Sgr may make 
some minor contribution to the VOD it is not the domi- 
nant origin of the stellar overdensity. 

8.4. SEKBO RRL 

An additional sample of RRL stars from the SEKBO 
(Moody et al. 2003) survey has recently been studied by 
Keller et al. (2008), with radial velocities presented by 
Prior et al. (2009). The Prior et al. (2009) subsample 
focuses on two groups of RRL which they group into the 
"VSS region" at a = 12.4-14 hr (i.e., A « 260° -290°) 
and the "Sgr region" at a = 20-21.5 hr (i.e., A 
10° — 40°). These two groups of stars are plotted against 
our ./V-body model in Figure [TO] (open triangles). 

Considering first the angular coordinate plot A Q vs. 
Bq, we note that neither of these groups of stars are 
well-centered on the predicted path of the iV-body tidal 
stream, lying at the extreme outer edge of the L1/L2 
streams and fully outside the angular path of the T1/T2 
streams. This is consistent with Figure 15 of Keller et al. 
(2008): The plane of the SEKBO survey intersects the 
Sgr plane at an angle and best overlaps at a ~ 18h (i.e., 
behind the Galactic Center) and a « 5h. In the regions 
traced by the Prior et al. (2009) subsample we confirm 
that the RRL tend to lie ~ 10° — 15° away from the bulk 
of the Sgr stream. It is perhaps unsurprising then that 
the distances (typically 16-21 kpc) and radial velocities 
of the RRL do not obviously match any particular branch 
of the Sgr stream. While Prior et al. (2009) discussed a 
possible match between some RRL and trailing Sgr tidal 
debris in the oblate-halo model of LJM05, no such match 
obviously exists for the updated model presented here. 

We note that the distance (~ 20 kpc) and radial ve- 
locity range (ugsr ~ —200 - +200 km s^ 1 ) of the 
A Q sa 10° - 40° sample of Prior et al. (2009) is similar 
to that found for a sample of red clump stars by Majew- 
ski et al. (1999; crosses in Figure IT9"|) . These red clump 
stars lie near the angular position of the LI stream, and 
LJM05 discussed the possibility that their large spread 



in radial velocities may indicate the range of velocities 
occupied by the Sgr stream at this longitude. Given the 
similarities in distance and radial velocity of the Prior 
et al. (2009) and Majewski et al. (1999) samples, it 
is at least likely that they trace the same Galactic sub- 
structure. If this substructure is confirmed to be the Sgr 
stream, it will suggest that the LI stream at A Q ~ 25° 
has a greater angular width, smaller distance, and larger 
spread in radial velocities than predicted by the current 
TV-body model. 

8.5. Subgiant Populations in the Leading Tidal Stream 

Keller (2009) has recently described a population of 
subgiant stars in the range a rj 120° — 180° from the 
SDSS that may correspond to the LI stream of Sgr. 
While this paper gives insufficient information to include 
these stars in Figure [TH we note that the A Q vs dis- 
tance trend illustrated in Figure 7 of Keller (2009) is 
reproduced by the revised ./V-body model to within ~ 1 
kpc, in contrast to previous V-body models of the Sgr 
stream (e.g., LJM05). Further characterization of this 
subgiant population may therefore provide even greater 
insight into the nature of the Sgr leading tidal stream. 

9. SUMMARY 

We have presented the first numerical model of the 
Sgr — Milky Way system that is capable of simulta- 
neously satisfying the majority of major constraints on 
the angular positions, distances, and radial velocities of 
the dynamically young tidal debris streams. In partic- 
ular, we have resolved the so-called "halo conundrum", 
whereby models were previously incapable of fitting both 
the angular position and distance/radial velocity trends 
of leading tidal debris. The publically-available JV-body 
model also makes concrete predictions for the location 
and kinematics of older tidal debris not yet conclusively 
observed in the Galactic halo. Future observational evi- 
dence for the presence/absence of such streams will place 
strong constraints on the total interaction time of the 
dSph with the Milky Way. 

The key ingredient in the success of this model is our 
introduction of a non-axisymmetric component to the 
Galactic gravitational potential, in the form of a nearly 
oblate dark matter halo with axial ratios 1.00/1.36/1.38 
whose minor axis lies within the plane of the Galactic 
disk plane within 7° of the line connecting the Sun and 
the Galactic center. While this inclined-halo explanation 
permits the model to reproduce the observational data, 
it is not strongly motivated within the current CDM 
paradigm, and may simply serve as a numerical crutch 
that mimics the effect of some as-yet unidentified alter- 
native (such as Sgr orbital evolution or the gravitational 
influence of the Magellanic Clouds) . 

Our TV-body model suggests that Sgr has been orbit- 
ing in the gravitational potential of the Milky Way for 
at least 3 Gyr and consists of a remaining bound mass 
Ms gr = 2.5t \'l x 1O 8 M . The total space velocity of Sgr 
is predicted to be t>sgr = 304 km s -1 , which at a distance 
-Dsgr = 28 kpc corresponds to an apparent proper motion 
of (fi a cosS, fig) — (—2.45,-1.30) mas yr _1 and a radial 
velocity of 142 km s _1 in the heliocentric rest frame. In 
contrast to some previous models, the Sgr stream is not 
expected to approach the solar position to closer than 
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~ 13 kpc. We have additionally demonstrated that it 
is possible to understand the observed evolution of the 
MDF along the Sgr tidal streams in terms of the dy- 
namical superposition of stars formed in multiple bursts 
of star formation and torn from the satellite at different 
times, suggesting that the chemical abundances of stars 
in Sgr at the present day may be significantly different 
than the abundances of those contributed to the Galactic 
stellar halo. 

Despite our success in fitting the strongest observa- 
tional constraints on the primary wraps of the leading 
and trailing tidal arms, it is uncertain what gives rise to 
the low surface brightness, high-declination 'bifurcation' 
of the stream that parallels the leading arm at ~ 15° 
higher declination in the SDSS footprint. We postu- 
late that substructure or anisotropic internal dynamics 
within the progenitor satellite may cause the observed 
bifurcation, indicating that future numerical models of 
the Sgr stream must explore the influence of the internal 
dynamics of the progenitor on the subsequently formed 



tidal tails. Such models must include realistic treatments 
of the dark and baryonic material as separate compo- 
nents and explore both rotation and orbital anisotropy 
of the visible matter. Moreover, the effects of evolution 
in the Sgr orbit and the peturbative gravitational influ- 
ence of Milky Way substructures such as the LMC must 
be investigated in greater detail. 
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